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Motivated by recent work of Lin, Balents, and Fisher [jjj, we compute correlation functions at zero 
temperature for weakly coupled two-leg Hubbard ladders and (N, N) armchair carbon nanotubes. In 
jlj it was argued that such systems renormalize towards the SO (8) Gross-Neveu model, an integrable 
theory. We exploit this integrability to perform the computation at the 5*0(8) invariant point. Any 
terms breaking the 5*0(8) symmetry can be treated systematically in perturbation theory, leading to 
a model with same qualitative features as the integrable theory. Using said correlators, we determine 
the optical conductivity, the single-particle spectral function, and the I-V curve for tunneling into 
the system from an external metallic lead. The frequency, w, dependent optical conductivity is 
determined exactly for u> < 3m (m being the fermion particle mass in the 50(8) Gross-Neveu 
model). It is characterized by a sharp "exciton" peak at ui = \f%m , followed by the onset of the 
particle-hole continuum beginning at u = 2m. Interactions modify this onset to a(uj + 2m) ~ uj 1 ^ 2 
and not the uj^ 1 ^ 2 one would expect from the van-Hove singularity in the density of states. Similarly, 
we obtain the exact single particle spectral function for energies less than 3m. The latter possesses a 
delta function peak arising from single particle excitations, together with a two-particle continuum 
for ui > 2m. The final quantity we compute is the tunneling I-V curve to lowest non- vanishing 
order in the tunneling matrix elements. For this quantity, we present exact results for voltages, 
V < (1 + v3)m. The resulting differential conductance is marked by a finite jump at lo — 2m, the 
energy of the onset of tunneling into the continuum of two particle states. Through integrability, 
we are able to characterize this jump exactly. 

All calculations are done through form-factor expansions of correlation functions. These give exact 
closed form expressions for spectral functions because the 5*0(8) Gross-Neveu model is massive: each 
term in the expansion has an energy threshold below which it does not contribute. Thus, we obtain 
exact results below certain thresholds by computing a finite number of terms in this series. Previous 
to this paper, the only computed form-factor of SO (8) Gross-Neveu was the two particle form factor 
of an 5*0(8) current with two fundamental fermions. In this paper we compute the set of all one 
and two particle form factors for all relevant fields, the currents as well as the kinks and fermions. 



I. INTRODUCTION 

Strongly interacting systems are a focal point of mod- 
ern solid state physics. The behaviour of these systems is 
often not captured by approximations, based on various 
'free particle' non-interacting models, since in the pres- 
ence of generically strong interactions the physics is typ- 
ically much richer, exhibiting qualitatively new features. 
Calculational tools to analyze the behaviour of strongly 
interacting theories are few in number. Amongst these 
are the tools of integrability, which are available for a 
number of special 1 + 1-dimensional models. These meth- 
ods allow for exact and detailed computations of many 
physically relevant properties of these theories. But as 
integrable models have to be of a special form, i.e. need 
to be 'fine-tuned', it is often believed that these tools 
are of little relevance for systems observed in the labo- 
ratory, since the latter would never generically be of the 
particular form necessary for integrability. 

However this need not be the case. Recently, in the 
context of transport experiments on tunneling point con- 
tacts in fractional quantum Hall effect devices H, it has 



been pointed out by |3| that when renormalization group 
(RG) thinking is combined with integrability, a special 
and 'fine-tuned' integrable Hamiltonian may actually de- 
scribe the behaviour of a realistic system. The Hamil- 
tonian of a generic system may be attracted under RG 
transformations to an integrable one, which then can be 
analyzed using the methods of integrability. In this way 
these powerful methods can be brought to bear on exper- 
imentally realizable systems and phenomena, and their 
predictions be directly compared with experimental data. 

In this paper we address another set of examples of a 
similar spirit. These are quasi one-dimensional interact- 
ing electronic systems, of which two prominent experi- 
mental realizations are two-leg Hubbard ladders (see for 
example |Q) and single- walled carbon nanotubes ||. 

Two-leg ladders have recently been the focus of much 
theoretical and experimental activity. At half-filling they 
are Mott insulators, exhibiting gaps to all excitations, 
and in particular a spin gap. These are typical examples 
of 'spin- liquids'. Upon doping, the gaps to all excitations 
except for those with charge-two survive [pj . The gapless 
charge modes induce quasi long-range superconducting 
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pairing correlations, with approximate <i-wave symmetry, 
reminiscent of underdoped cuprate superconductors. 

Carbon nanotubes are novel materials whose mechan- 
ical and electronic properties promise potential for new 
technological applications || . They are formed by wrap- 
ping graphite sheets into cylinders of nanoscale dimen- 
sions. They support electronic excitations, which, for 
a prominent member of the nanotube family, the arm- 
chair (n, n) type, can be described by the same theo- 
retical model as that used for the two-leg Hubbard lad- 
ders [^J|]. Even though these systems would be one- 
dimensional band metals in the absence of interactions, 
they become Mott insulators at half-filling due to the 
presence of short-ranged electronic interactions, which 
play an important role due to their one-dimensional na- 
ture. It is these interaction effects that we analyze ex- 
actly in this paper using the powerful methods of integra- 
bility. After experimental techniques had been developed 
to fabricate long single- walled nanotubes with high yields 
in the laboratory, this field of material science has seen 
a explosive development ||. Electronic properties can 
be measured relatively easily by attaching metallic leads 
or by tunneling into these materials with scanning 
tunneling microscope (STM) tips. The practical feasibil- 
ity of such tunneling experiments from an STM tip into 
an individual single- walled nanotube placed upon a gold 
substrate (screening long-range Coulomb forces) has re- 
cently been demonstrated by J.W.G. Wildoer et al. 
Thus, the predictions that we make in this paper, es- 
pecially those for the tunneling /(V)-curve, have direct 
bearing on possible future experiments on these materi- 
als. 

As said, both systems, two-leg ladders and single- 
walled armchair nanotubes, would be one-dimensional 
metals in the absence of electron interactions. These 
are described theoretically (on scales much smaller than 
the non-interacting band width) by two species of spin- 
ful massless Dirac fermions in (1 + 1) dimensions 1 . For 
the ladder compounds, the two species arise in an obvi- 
ous way from the two rungs of the ladder, whereas for 
the nanotubes they arise from the particular band struc- 
ture of the underlying hexagonal graphite lattice, char- 
acterized by two Fermi points in the Brillouin zone |^,^| . 
These massless Fermi surface excitations interact with 
short-range interactions whose detailed nature is deter- 
mined by non-universal microscopic considerations. 

A notable observation was made recently by Lin, Ba- 
lents, and Fisher jlj]. These authors argued that within 
an 1-loop RG any such model with generic, non-chiral, 
short range interactions flows at half-filling into a the- 
ory with an immense symmetry, namely the 50(8) sym- 



1 These two species have equal Fermi velocities due to 
particle-hole symmetry present at half-filling. 



metric Gross-Neveu model. This model not only has a 
large 50(8) global symmetry, which encompasses an one- 
dimensional version of SO (5) recently advocated by S. C. 
Zhang [ pd| , but in addition has an infinite number of hid- 
den conservation laws, which are a consequence of the in- 
tegrability of this model. The 5*0(5) subgroup symmetry 
is the same as that studied by Jl2| in the context of lattice 
models of interacting two-leg ladders. We point out that 
it would be straightforward to construct a corresponding 
50(8) invariant lattice model of a two-leg ladder. 

The use of the RG in § can be understood in the fol- 
lowing sense. Since the analysis in JO is based on an 1- 
loop RG, the initial microscopic (bare) interactions must 
be small enough so that the integrable 50(8) invariant 
RG trajectory is approached sufficiently closely after a 
number of RG steps, before leaving the range of validity 
of the 1-loop RG equations. Whenever this is the case, 
it is argued that the integrable model is approached in- 
dependently of the (sufficiently weak) values of the bare 
interactions. The situation for the 2-leg ladder is thus 
similar in spirit to that of the point contact device en- 
countered in H), where only a single operator was rele- 
vant, and this relevant operator was integrable. In the 
latter case all other interactions were irrelevant in the RG 
sense, and could in principle be treated perturbatively. 

The requirement of the RG that the interactions be 
short-ranged is natural in the case of the Hubbard lad- 
ders. However it may not seem so in the case of the 
carbon nanotubes. Recent theory |l3| and experiment 
have discussed the case where long-ranged Coulomb 
forces drive Luttinger liquid behaviour in single-walled 
carbon nanotubes. However we do not have such situ- 
ations in mind for the paper at hand. Rather we want 
to consider situations such as those found in the experi- 
ments of J.W.G. Wildoer et al. ||] where the long range 
forces are screened. 

Although the restriction to such experiments in the 
case of the carbon nanotubes places us upon safe ground, 
it is not inconceivable that experiments where the long- 
ranged forces are present would nevertheless see be- 
haviour indicative of the SO (8) symmetry. An un- 
screened force translates into an unusually large bare 
coupling (in comparison with other bare couplings) in 
the forward scattering direction. However this does not 
mean the RG is inapplicable. The RG still indicates a 
potential enhancement in the symmetry. Because of the 
large bare coupling, the RG must be run a longer time be- 
fore any enhancement would be seen but nevertheless an 
enhancement may well occur at some low energy scale. In 
terms of the experiments in |j"4[ , this would mean that at 
medium energy scales, Luttinger liquid behaviour would 
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predominate, while at much lower energy scales, SO (8) 
behaviour would be expected. However at current stand- 
ing, the material science is not advanced to the point 
where it is possible to accurately probe the very low en- 
ergy behaviour. But the potential for advancement in 
this area is ever present. 

The RG analysis further requires the bare couplings 
to be weak. With Hubbard ladder compounds, this con- 
dition will not be generically met, although it certainly 
will not be universally violated. However with (N,N)- 
armchair carbon nanotubes, the bare couplings are natu- 
rally weak. It is one of the hallmarks of the physics of the 
(N,N) armchair carbon nanotubes that the electrons are 
delocalized around the circumference of the tube. This 
in turn leads to a scaling of the effective short-ranged 
interaction by 1/N, making it naturally small Q|. 

It can, however, be questioned on a more fundamen- 
tal level whether an 1-loop RG adequately describes the 
system's behaviour. Difficulties with the analysis in fj] 
take two forms. As a first objection, the authors of 
point out that an RG flow can imply a symmetry restora- 
tion which in fact does not occur. As an example they 
consider a U(l) symmetric Thirring model, 

C = + \ 9] \{j z ) 2 + \gx[{k? + (jy)% (LI) 

where j£ — ^^^a a ^. Although the 1-loop RG equa- 
tions for this model seems to indicate a generic sym- 
metry restoration to a more symmetric SU(2) case (i.e. 
g± = <7y ), this in fact only occurs in a certain region of 
coupling space. For n — \g±\ > —gn > \g±\ > 0, the U(l) 
model maps onto the sine-Gordon model with interaction 
cos(f3(f>) fill , where (3 is given by 

(3 2 = 8-k — 8/x; [i — cos -1 [cos(g||)/ cos(gj_)]. (1.2) 

The value of (3 completely characterizes the model. While 
g_i_ and g\\ flow under the RG, the particular combination 
of these parameters forming /3 does not. Thus for this 
particular region of parameter space the model moves no 
closer to the SU(2) symmetric point under an RG flow. 
In other regions however (for example \g±\ > \g\\\), the 
situation is better; the effect of the anisotropy in the cou- 
plings is exponentially suppressed. 



However it is reasonably clear that such pessimism is 
not warranted in the analysis of the RG of 0]. A salient 
criticism of |ll| is that in considering the action of the 
renormalization group, they fail to consider the conse- 
quences of working in the scaling limit. The scaling limit 
is exactly the limit in which a field theory becomes avail- 
able. In turn, the scaling limit places constraints upon 
the possible range of bare couplings consistent with a 
field theory. In the case of sine-Gordon, the underlying 
integrability/solvability of the theory allows explicit in- 
vestigation of this question. It is found that the allowed 
range is such that even moderate anisotropic deviations 
are forbidden The scaling limit, in other words, en- 
forces isotropy. In cases where there are RG flows indi- 
cating an enhancement in symmetry, this turns out to be 
a general phenomena and it leads to an expanded notion 
of symmetry restoration fl7|j . 

On a more concrete level, the breaking of the SU(2) 
symmetry considered in []l5f is a rather special case. Ul- 
timately, the parameter p in the sine-Gordon model is 
protected under an RG flow by the presence of a quan- 
tum group symmetry arrived at by deforming a Yan- 
gian symmetry present at the SU(2) point. There is, 
however, no such known way to deform the Yangian in 
50(8) Gross-Neveu. Indeed the natural generalization of 
the sine-Gordon model to SO(8) is not to SO(8) Gross- 
Neveu but to an affine toda SO(8) theory where such a 
deformation of the Yangian symmetry is possible jl8| . 

Another question that one must ask in looking at the 
analysis in |l5f| is how the choice of the symmetry break- 
ing terms affects the symmetry restoration. The sine- 
Gordon model still possesses a U(l) symmetry. How- 
ever it is certainly possible to consider perturbations that 
break this U(l). Such perturbations would destroy the 
quantum group symmetry of sine-Gordon and thus might 
lead to symmetry restoration. This would be perhaps 
closer to the RG analysis of (!]] where a large number 
(nine) of marginal perturbations were included. We in 
fact consider exactly such a situation in |17j and find 
that indeed there is symmetry restoration. 

The second objection to the analysis of |lj is its omis- 
sion of chiral interactions that alter the Fermi velocities 



2 We note in passing that the authors of |15| consider an 
anisotropic Gross-Neveu model, a model of direct relevance 
to the situation at hand. They conclude through a mean- 
field/large N limit computation that the model is intrinsi- 
cally anisotropic thus throwing doubts upon the analysis in 
Lin et al. jl|. There appear to be several severe problems, 
however, with their analysis. The most telling is the manner 
in which the bare couplings are scaled in the large N limit. 
The anisotropic model they consider has three bare couplings. 
One is chosen to not scale at all, one scales as 1/N, and the 
last scales as 1/N d/2 , d < 1. With this scaling (Ta the model 
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p9| . Such interactions, although they are absent from 
the 1-loop RG, likely play a role at higher order. How- 
ever their effect is less drastic than envisioned in [ p^[ . 
There a scenario was considered where the invariant RG 
trajectory of higher symmetry was inherently unstable to 
perturbations. However the SO (8) RG ray in Q has a 
basin of attraction of finite measure. The effect of chi- 
ral interactions is to then slightly alter the direction of 
the ray [^0). In turn the ratio of masses of the various 
excitations are slightly perturbed away from one. 

In taking account of these objections, prudence sug- 
gests a modification in the understanding of the RG anal- 
ysis of {|] . This analysis in Hj and ^cj tells us that while 
the RG flow does not restore an exact symmetry, it leaves 
us close to the symmetrical situation. In particular, it in- 
dicates that while the masses in the actual system may 
differ from their 5*0(8) values, they do not wildly di- 
verge. One then understands the 50(8) Gross-Neveu 
theory, not as precisely representative of the actual sys- 
tem, but in near perturbative vicinity of it, that is, as an 
excellent starting point about which to perform pertur- 
bation theory in the non-integrable interactions breaking 
50(8) in much the same spirit as done for a non-critical 
Ising model in the presence of a magnetic field plj . 

In doing perturbation theory about 50(8) Gross- 
Neveu it is important to emphasize two salient points. 
Firstly 50(8) Gross-Neveu already contains all the ba- 
sic features of the strongly interacting electron system. 
It thus provides a much better starting point for pertur- 
bation theory than the two band metal from which it 
arose. In contrast, perturbation theory about the two 
band metal could not hope to capture, even qualita- 
tively, characteristics of the interacting system. Secondly, 
50(8) Gross-Neveu is a massive theory and so poses none 
of the attendant problems of perturbation theory with a 
massless model (as a two band metal is). The pertur- 
bative series is well controlled; the perturbed theory is 
connected in a continuous fashion to 50(8) Gross-Neveu. 

In this paper we focus on the integrable model at half- 
filling and at zero temperature 3 and as such this work 
should be considered as the starting point for the analy- 
sis of a perturbed 50(8) Gross-Neveu ladder /nanotube. 
(We do, however, in Section 3 make a rough sketch of the 
effect of perturbations on our results. A more detailed 
analysis awaits further work.) Some consequences of the 
integrability of the theory have been exploited in p) , no- 



tably those following directly from the structure of the 
elementary excitations ('particles') of the theory. If U is 
the typical strength of the interactions and t a measure of 
the bandwidth of the non-interacting model, the Gross- 
Neveu interactions generate dynamically a single mass 
scale m oc te^ l ^ u . Due to the integrability, the masses 
of all excitations are universally related to m, universal 
mass ratios being known exactly. This, together with the 
knowledge of the quantum number assignment to these 
excitations, has been exploited in [Q to make, inter alia, 
a number of predictions for correlations functions, largely 
of a qualitative nature. 

Due to its integrability, an immense amount of infor- 
mation is known about this strongly interacting system. 
In this paper, we put this information to use, to com- 
pute the exact functional form of spectral functions, at 
zero temperature and on scales below certain thresholds 
(to be discussed below). To this end we use the form 
factor approach described in [p2| . In particular, we have 
obtained exact expressions for the optical conductivity, 
the single-electron spectral function (experimentally ac- 
cessible via photoemission experiments), as well as the 
tunneling density of states (experimentally observable for 
example by measurement of the differential conductance 
for tunneling from a metallic lead into a nanotube (^]). 

We emphasize that our results are exact on energy 
scales below certain thresholds. In order to appreciate 
the significance of this, it is important to recall that even 
in an integrable system, correlation functions cannot in 
general be computed exactly. However if a mass gap is 
present for all excitations (such as for the Gross-Neveu 
model) exact results for spectral functions can be ob- 
tained below certain thresholds. This can be easily un- 
derstood from the basic ideas underlying the form factor 
approach to correlation functions in integrable systems, 
which we now briefly review. The key feature of an inte- 
grable system is the exact knowledge of a basis of eigen- 
states of the fully interacting Hamiltonian. At the root 
of integrability is a well defined notion of "particles" , or 
"elementary excitations" in the fully interacting system. 
These particles scatter off each other only with two-body 
5-matrices, that is, all particle production processes are 
absent and particle number is conserved. This is due 
to special conservation laws which exist in an integrable 
model, preventing the decay of these particles. In this 
sense, an integrable system is similar to a Fermi liquid. 



is anisotropic. However all that was done is to examine a 
system with a diverging bare anisotropy and conclude that 
the effective theory is similarly anisotropic. And nowhere is 
the RG, which in this model indicates symmetry restoration, 
allowed to act. As such, we believe this example has little 
bearing on the situation at hand. 
3 The doped system is examined in M . 
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An additional feature is that new particles can arise as 
bound states of already existing ones. However the to- 
tal number of different types of particles is finite which 
makes the system analytically tractable. For the 50(8) 
Gross-Neveu model these particles consist of fermionic 
particle states of mass m, namely one octet of funda- 
mental fermions and two octets of kinks. The excitations 
of the ladder or nanotube, possessing the quantum num- 
bers of the electron, are represented by eight of these 
kinks. In addition there are 29 bosonic particle states 
organized as a rank-2 SO(8) anti-symmetric tensor and 
a singlet, all of mass y/3m. These are bound states of 
the fermionic mass-m particles. The exact eigenstates of 
the interacting Hamiltonian arc simply n-particle states, 
\n;s n ), where s n collectively denotes the momenta and 
species labels of the n particles 4 . 

The form-factor representation of any correlation 
function is obtained by inserting a resolution of 
the identity corresponding to the basis of parti- 
cle eigenstates. Thus, for an operator, 0(x, r), 
we write the spectral decomposition schematically (r 
denotes imaginary time, and T time ordering) 5 



G%(x,t) = (0\T(O(x,t)O^Q,0))\0) 

OO 

^TJ2(0\O(x,0)\n;s n )e- rE >»x 

(n; S „|0t( O ,O)|O), (r>0), (1.3) 

where E Sn is the energy of the eigenstate, |n;s n ). In 
an integrable model the matrix elements of the physical 
operator between the vacuum and the exact eigenstates 
can in principle be computed exactly from the 2-body 
S-matrix. However the calculation of these matrix el- 
ements, as well as the evaluation of the sums/integrals 
^2 S , becomes increasingly cumbersome as the particle 
number n becomes large, so that the full expression for 
the correlation function cannot be evaluated in closed 
form. Often, however, a truncation of the sum at the level 
of two or three particle states already provides a good ap- 
proximation to the full correlation function |23| On 
the other hand, this truncation is no longer necessary 
in a massive theory, if one considers the corresponding 
spectral function. Only eigenstates with a fixed energy, 
u>, contribute to the spectral function: 
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1 °° f 
ImG?(x, -iu + 8) = \ (0\O(x,0)\n; s n ){n; s n \O\0, 0)|0)J(w - E S J 

71=0 S„ ^ 

-e<O|0t(0, 0)|n; s„)(n; s n \0(x, 0)\0)8(u + E Sn )\, (1.4) 



where e = ± for fields, O, that are bosonic/fermionic. 

Since in a massive theory the creation of an extra par- 
ticle in the intermediate exa ct ei genstate costs a finite 
amount of energy, the sum in 1.4 is finite. For example, 
when u> is smaller than the energy of all three-particle 
states (i.e. when uj is below the three-particle thresh- 
old), then only the form factors with one and two par- 
ticles (n = 1,2) have to be determined. This is what 
we have done in this paper for all the spectral functions 
we computed. Thus, the results obtained from the form 
factor method for spectral functions are exact in massive 
theories. 

The results for the physical quantities we compute, 
namely the optical conductivity, the single particle spec- 



tral function, and the tunneling I —V curve, are summa- 
rized in Section 3. The basic features of these results are 
as follows. The optical conductivity, Re[<r(a;)] (so called 
because it would be measured in a reflectivity experi- 
ment) has two notable features: a delta function peak at 
ui = \/3m corresponding to an excitation of one of the 28 
rank two tensorial bosonic particles, and a continuum of 
two particle states beyond u = 2m (see Figures 1 and 2). 
In a free theory, the van-Hove square-root singularity in 
the density of states at the two particle threshold would 
lead Re[cr(u;)] to behave as 1/y/ui — 2m 6 . However the 
current matrix elements vanish at threshold changing the 
behaviour to \/lu — 2m. This vanishing is ultimately the 
result of generic interactions becoming strongly renor- 



4 This can be made more explicit by writing, |n; s n ) — 
|(pi,oi); (p2,a 2 ); (p„,a„)). 

5 The sum over s n is meant to include integrals over the mo- 
menta of all particles. 

6 This is the result found in jjj. In this work the optical 
conductivity is computed in the large TV limit of SO(2N), in 
effect an RPA approximation. In this limit the theory consists 
of four massive but non-interacting Dirac fermions. 
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malized at low energies |2!| . However its exact behaviour 
as one moves away from threshold can only be extracted 
through integrability. 

As the SO (8) Gross-Neveu model characterizes the low 
energy behaviour of the Hubbard ladders/carbon nan- 
otubes, we are able to describe the single particle spec- 
tral function only near the Fermi points. Here we find 
two principal features. There is a sharp peak (infinitely 
so at T = 0) describing the single electron contribution 
to the spectral function. As the system is interacting we 
expect to find a continuum of higher particle contribu- 
tions to the spectral function. In physical terms these 
higher particle contributions would take the form of neu- 
tral particle-hole excitations together with an electron. 
In the Gross-Neveu language said excitations take the 
form of a kink (akin to the electron) and a fundamen- 
tal Gross-Neveu fermion (akin to the particle hole exci- 
tation). The threshold of this two-particle contribution 
occurs at an energy, lo = 2m. 

One of the quantities most easily measured in experi- 
ments on nanotubes is the tunneling density of states. 
This quantity can be measured by tunneling from a 
metallic lead, such as an STM tip, into the nanotube. 
When applying a voltage bias, a current will flow. This 
current is a non-linear function of the applied bias. To 
lowest non-vanishing order in the tunneling matrix ele- 
ment, this current can be related to the single-particle 
spectral function at energy lo = V (using standard rea- 
soning). This relationship is non-perturbative in the 
voltage V. Differentiating the so-obtained non-linear 
current-voltage characteristic, the differential conduc- 
tance for tunneling into the nanotube is easily seen to 
equal the spectral function of the single-particle Greens 
functions. Physically, it is a 'spectroscopic probe', or 
'measure' of the density of states in the fully interacting 
nanotube at energy V, even though it is not equal to it as 
the spectral function is equal to the density of states mul- 
tiplied by the form-factor matrix elements. A look at Fig- 
ure 5 reveals an interesting novel feature, entirely due to 
the interactions, and not visible in a band-semiconductor: 
the differential tunneling conductance develops a jump 
of a finite magnitude at the two-particle threshold. This 
might have been expected on physical grounds, since at 
this energy extra states suddenly become available for 
current transport. Nevertheless, is it by no means ob- 
vious that this would correspond to a discontinuity in 
the tunneling density of states; the two-particle channel 
could also open up gradually. This depends on the de- 
tailed behaviour of the matrix elements near the thresh- 
old, which we have computed exactly from integrability. 
Our analysis shows that these matrix elements behave in 



such a way that there is a finite jump. 

The one and two particle form factors have not in gen- 
eral been computed before for the SO (8) Gross-Neveu 
model. The sole exception is the matrix elements of the 
current operators with the fundamental fermions [j30| . 
The calculation of the remaining form factors, and in 
particular those involving the kinks, makes up the tech- 
nical part of this paper. To compute the form factors, we 
employ a series of algebraic constraints that arise from 
consistency with known two body scattering matrices, 
unitarity, Lorentz invariance, and braiding relations. The 
first three are common in form factor calculations (see for 
example p6| p8 23 2^]3(|). The latter is perhaps more 
unusual [ 27| , ^2|j3j| ^ Braiding relations arise as general- 
ized commutators of fields, ip a : 

r(x, W b (y, t) = R^ d {y, t)r(x, t), x > y, (1.5) 

where is termed the braiding matrix. They are both 
indicative of the non-locality of the fields for which we 
compute form factors and of exotic symmetries. Non- 
trivial braiding relations are tied to the existence of quan- 
tum group symmetries |32] ]. In this case, the relevant 
quantum group symmetry is an SO (8) Yangian. This is 
similar to the situation found in the well-studied sine- 
Gordon model at its SU(2) point [Q, where an SL(2) 
Yangian symmetry is present. 

The SO(8) Yangian is the operative symmetry of the 
model. The particles in the theory are thus organized 
in terms of its finite dimensional representations and not 
SO (8) 's. There is however a large degree of correspon- 
dence between the two sets of representations. The three 
fundamental eight dimensional representations of SO (8) 
are also irreducible representations of the Yangian. How- 
ever the 28 dimensional second rank anti-symmetric ten- 
sor representation does not appear as a representation of 
the Yangian. Rather under the Yangian it is combined 
with the one-dimensional scalar representation into a sin- 
gle 29 dimensional representation 7 . Thus, unsurprisingly 
the scalar particle together with the 28 particles of the 
anti-symmetric tensor share the same mass, \/3m. 

The paper is organized as follows. In Section 2, we re- 
view the arguments given by Lin, Balents and Fisher pj 
showing how the various massive phases of the Hubbard 
ladders/nanotubes are related to an 50(8) Gross-Neveu 
model. In doing so we establish notation. In Section 
3, as stated previously, we summarize our results for the 
optical conductivity, the single-particle spectral function, 
and the tunneling I-V curve. Sections 4 and 5 are devoted 
to computing the form factors. For readers uninterested 
in the details of this analysis, the last part of Section 5 



We would like to thank N. MacKay for stressing this fact 
to us. 
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provides a summary of results. Specifically, Section 4 re- 
views the S-matrices and attendant group theory for the 
50(8) Gross Neveu model, while Section 5 presents the 
actual form factor computations. 



II. HUBBARD LADDERS TO SO (8) 
GROSS-NEVEU 

Here we will briefly review the connection between 
Hubbard ladders (and related armchair carbon nan- 
otubes) and the 50(8) Gross-Neveu model developed in 
jjj . Specifically, we summarize the map between the two 
models, and interpret the excitation spectrum and fields 
of the 50(8) Gross-Neveu model in terms of the original 
ladder model. 



A. Map: D-Mott Hubbard Ladders to SO {8) 
Gross-Neveu 

The weakly interacting Hubbard ladder has five differ- 
ent phases, one massless and four massive §. The four 
massive phases, D-Mott (a spin liquid with approximate 
short-range ei-wave pairing symmetry), S-Mott (a spin 
liquid with approximate short-range s-wave pairing sym- 
metry), spin-Peierls (SP) (electrons are dimerized along 
the legs of the ladder), and a phase with charge density 
wave order (CDW), correspond to various combinations 
of attractive and repulsive interactions. We first focus 
upon the D-Mott phase which is characterized by generi- 
cally repulsive interactions. Once we discuss the D-Mott 
phase in detail, we will return to the other three massive 
phases and demonstrate that each has a distinct 50(8) 
symmetry In the massless phase, termed C2S2 (i.e. 
two charge bosons and two spin bosons), interactions are 
irrelevant, and so the phase has two trivial independent 
50(8) chiral symmetries. For this reason we will not 
consider it here. 

We follow H in relating the D-Mott Hubbard ladder 
to an 50(8) Gross-Neveu model. We begin with non- 
interacting electrons hopping on a ladder: 



H = -}Xta\ a (x + l)a la (x) + ta\ a (x + l)a 2a (x) 

+ t ± a{ a (x)a2 a (x)+h.c). (2.1) 

Here the ai/aj are the electron annihilation/creation op- 
erators for the electrons on rung I of the ladder, x is a dis- 
crete coordinate along the ladder, and a —1,1 describes 
electron spin, t and ij_ describe respectively hopping be- 
tween and along the ladder's rung. 

The first step in the map is to reexpress the a's of Hq 
in terms of bonding/anti-bonding pairs: 



With this transformation, the Hamiltonian can be di- 
agonalized in momentum space in terms of two bands. 
Working at half filling, particle-hole symmetry dictates 
that the Fermi velocities, v F j, of the two bands, j = 1, 2, 
are equal. As we are interested in the low energy be- 
haviour of the theory, the Cj a 's are linearized about the 
Fermi surface, kpj\ 



(2.3) 



where L, R corresponding to the right and left moving 
modes about the Fermi surface. With this Hq becomes, 

H = v F J dx^2 [cR ja idxCR ja - c\ ja id x c L ja] ■ (2.4) 

ja 

The next step is to bosonize the c's: 

c P]a = K ja e 1 ^ , P = +, - = R, L. (2.5) 
Here Klein factors satisfying 

{Kj a , Hif3} = 2SijS a f3. (2.6) 

In terms of these four Bose fields, four new Bose fields 
are defined (effectively separating charge and spin): 



5P1 



>P2 = - 



?P3 



9p 4 



opii ■+■ q>p\i -t q>p2\ -t q>P2i) 

6pif — 0P1X + 4>P2\ — 4>P2l) 
6pi| — 4>P1{ — (j>P2\ + 4>P2l) 
0P1T + 4>Pll — 4>P2] — 0P2J.) 



(2.7) 

Note that <fip4 has a relative sign between the right and 
left movers. This sign effectively masks the 50(8) sym- 
metry o f th e original Hamiltonian. Without the sign 
change, 2.7 is no more than a triality transformation 
(see Section 2.4). If we then refermionize with these new 
bosons, i.e., 



1. 



Cj a = —7=(aia + (-l)- 7 a 2 a). 

V2 



(2.2) 



*Pa = K a e i *°, 

*P 4 = PK 4 e^ 4 , 
where the Klein factors are given by 

we find for the free Hamiltonian, 



(2.8) 



(2.9) 
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H= f dx^[*l L id x * aL -*i R id x y a R], (2.10) 

where the Fermi velocity, vf, has been set to 1. 

The point to the bosonization, change of basis, and 
refermionization only becomes apparent when one con- 
siders interactions. The discovery in |lj was that if one 
writes down a generic set of weak, left-right, repulsive in- 
teractions between electrons on the rungs, expresses these 
interactions in terms of the refermionizcd fcrmions, and 
then allows the couplings to flow under an RG, one finds 
that the interacting Hamiltonian is attracted to, 

4 

Hint = ff[y)(i*Li*afl " i*ifl*ai)] 2 - C 2 " 11 ) 



a=l 



This is, of course, H int for the 50(8) Gross-Neveu model. 

It will sometimes prove convenient to recast the the- 
ory in terms of Majorana fermions, ij) a p. In terms of the 
Dirac fermions, $„p, they are given by 

^ a p = ^=(ip 2a ,p + i^2a-i,p), (a=l,...,4). (2.12) 



In this basis, Hint can be recast as 

H int =gG a B b Gt, (a > b = !,..,<< 



(2-13) 



where G P b — iifj a pipbP is one of the 28 50(8) Gross- 
Neveu currents. 



B. D-Mott Excitations in Gross-Neveu 

The Gross-Neveu SO(8) model has an exceedingly rich 
spectrum. There are 24 fermionic particles of mass m or- 
ganized into one eight dimensional vector representation 
and two eight dimensional spinor representations. We 
denote the particles of the vector representation by A a , 
a = 1,...,8. The A a 's are the Majorana fermions of 
2.12. The kink particles, in turn, will be denoted by A a . 
Here a is of the form a = (±1/2, ±1/2, ±1/2, ±1/2) and 



so takes on 16 values. These 16 particles decompose into 
the two eight-dimensional spinor representations. This is 
discussed in more detail in Section 4. 

Beyond the eight dimensional representations, there 
are 29 bosonic particle states of mass -\/3™, transforming 
as a rank- two tensor of dimension 28 and a singlet. To- 
gether they form a representation of the 50(8) Yangian 
symmetry. These particles can be thought of as bound 
states of either two kinks or two fundamental fermions. 

As 50(8) is a rank 4 algebra, the 50(8) Gross-Neveu 
model has four Cartan bosons (i.e. the cf)p a , a = 1, . . . 4) 
and so its excitations are characterized by four quantum 
numbers, Ni, i = 1, . . . , 4. With the Majorana fermions, 
the combination 



carries quantum number N a = ±1, Nb = 0, b ^ a. The 
quantum numbers carried by the kinks A a are directly 
encoded in a. If a = (01,02,03,04), a,j = ±1/2, the A a 
carries the quantum numbers, Ni = a^. The quantum 
numbers carried by the rank two tensor states can be 
directly deduced from the particles forming the bound 
state. As we will always think of the bound states in this 
way, we will not list their quantum numbers directly. 

The last thing needed in the section is to identify 
the relationship between the quantum numbers, Ni, and 
the physical quantum numbers of the system, the z- 
component of spin, S z , the charge, Q, the difference 
in z-component of spin between the two bands, 5i2, 
and the "relative band chirality", P12, defined as P12 = 
Nri — Nli — Nr2 + Nl2, where Npj is the number elec- 
trons in band j with chirality P. In p| it was found: 



{Ni = 1,0,0,0) <-> 


(Q 


= 2,5, 


= 0, 5 12 


= 0, Pa 


= 0); 


(0,N 2 = 1,0,0) 


(Q 


= 0,5, 


= 1, 5i2 


= 0, P12 


= 0); 


(0,0,^3 = 1,0)^ 


(Q 


= 0,5 Z 


= 0, 5 12 


= l,Pl2 


= 0); 


(0, 0, 0, Ni = 1) <-> 


(Q 


= 0,5 2 


= 0, 5 12 


= 0, P 12 


= 2). 



(2.15) 

With this assignment, we can see that the vector repre- 
sentation of fundamental fermions corresponds to states 
of two electrons in the original formulation. For example, 
the fermion A2 ± %A\ carries charge ±2 and no spin (the 
cooperons), and the fermion A^±iAs carries spin, 5 2 = 1, 
and no charge (the magnons). This makes concrete the 
earlier comment that the change of basis of bosons in 2.7 
affects a charge-spin separation. The spinor represen- 
tations, the kinks, in turn correspond to single particle 
excitations as their quantum numbers are combinations 
of Ni/2. 



C. D-Mott Gross-Neveu Fields 

In this section we make contact between the fields of 
the 50(8) Gross-Neveu model and the original fields of 
the Hubbard ladders. As we have already discussed, the 
fundamental (Dirac) fermions of the vector representa- 
tion are given by 



* aP — ri a c , 



*aP = PK a e 



i<PaP 



(2-16) 



A 2a ± iA 



2a— 1) 



(2-14) 



and carry quantum numbers corresponding to two elec- 
tronic excitations. However the $„p are fermionic, 
whereas such excitations are bosonic. As such, ^ aP are 
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not simply related to a fermionic bilinear of the original 
electrons but must be a fermion bilinear multiplying some 
non-local field (a Jordan- Wigner string). As we will not 
compute correlators involving such fields in this paper, 
we will not elaborate upon this (see Q). 

As discussed previously, the kinks correspond to sin- 
gle particle excitations. Thus we expect to find that the 
kink fields are related to the original electron operators. 
This is true in part. There are 32 kinks in total (counting 
both left and right movers), but only sixteen electron op- 
erators, the c's and c^'s (four for each of the four Fermi 
points). So we expect only 1/2 of the kinks to correspond 
to actual electron operators. 

We represented the fundamental fermions in terms of 
the four Cartan bosons. There is a corresponding repre- 
sentation for the kink fields: 



Ja-4> r 



(2.17) 



where <j) — (0i, <f> 2 , <t>3 , <Pa ) ■ The kink fields that then cor- 
respond to the electron operators c's are as follows: 



D. Other Massive Phases 

By adjusting the signs of interaction couplings, three 
other massive phases, S-Mott, SP, and CDW, can be ob- 
tained jlj. Each of these phases is characterized by a 
distinct 50(8) symmetry. However they share a com- 
mon 50(5) subgroup. The 50(5) symmetry is the one- 
dimensional analog of the 50(5) symmetry recently pro- 
posed by S.C. Zhang as a means to unify antiferromag- 
netism and superconductivity. Each of these new 50(8) 
symmetries is readily expressible in terms of the D-Mott 
50(8) symmetry through considering the four Cartan 
bosons <&p a As the first two Cartan generators be- 
long to the 50(5) subgroup they, however, remain un- 
changed. It is the latter two bosons that are affected. 
For completeness we review each of new phases in turn. 

S-Mott: 

The defining four Cartan bosons of the S-Mott phase, 
4>p a , are related to those of the D-Mott phase via 



crit 

CR2-[ 
CR21 

cm i 



i(4>lR+4>2R+4>3R + <i>4R) /2 . 

) 

i((/>lR+<p2R — 4>3R-4>>lR)/2- 
1 

i(4>lR—<t>2R+'l>3R-4>4R)/2. 

) 

— 4>2R-4'3R + <t>4.R)/2. 



(even chirality) 



clit 

C-L2] 
C L2[ 
CLll 



i(4>lL+4>2L+<t>3L—<p4L)/'2. 

; 

«(01-L+02i-03t+04i)/2. 

5 

i(4>lL -<j>2L + 4'3L+<p4L)/2. 

) 

-<t>2L—<t>3L—<t>4.L)/2. 



(odd chirality). 



(2.18) 



With hcrmitian conjugates, this totals to sixteen fields. 
The ~ sign is meant to indicate that these equivalences 
hold up to Klein factors. The cpj a 's, of course, are 
fermionic. However the kink fields as defined are not. 

The last set of fields that are of concern to us are the 
currents. The electric current of the ladder has the lattice 
representation 



J i 



la L 



aj a {x)ai a (x + 1) - a\ a {x + l)ai a (x) 



(2.19) 

where we have summed over the contribution coming 
from each spin (a) and each leg (/) of the ladder. Taking 
the continuum limit, J equals, in Gross-Neveu language, 

J~ isinkpi d t (j)i ~ G 12 , (2.20) 

where G\ 2 is one of the 50(8) currents discussed in [2.13 



^P4 



1,2,3 



Ptt/2. 



(2.21) 



The sole difference between the 50(8) algebra of the S- 
Mott phase and that of the D-Mott phase is a shift of 
<pPa by Pn/2. This shift forces a sign change in the pair 
field correlator, changing the symmetry from d-wave to 
s-wave. The shift, however, does not change the exci- 
tation spectrum or its attendant assignment of quantum 
numbers, nor does it change field assignments beyond a 
phase. 

Spin-Peierls (SP): 



The defining four Cartan bosons of the SP phase, 
are related to those of the D-Mott phase via 



kSP 



kSP 
°Pa 



lSP 
3 P3 



1,2,4 



P<t>ps; 



(2.22) 



that is, the sole change is to flip the sign on The ef- 
fect of this sign change is two-fold. The quantum number 
associated with the third Cartan boson is now defined to 
be 

N 3 = 5i2 = -(JVbti - Nmi - N R]2 + N R[2 )- 

-(N Ln -N Lll -N Ln + N Ll2 ), (2.23) 



where Npj a is the number of electron of chirality P in 
band j with spin cr, and so A3 is the relative right- 
moving spin between the two bands minus the relative 
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left-moving spin between the two bands. Moreover the 
left moving electrons, cl, are now to be identified with 
kinks of even chirality as opposed to the odd chirality 
kinks for the D-Mott phase (see 2.18). 



Charge Density Wave (CDW): 



The defining four Cartan bosons of the CDW phase, 
t>Pa W , are given by 



• CDW 



lCDW 
9 P3 



lCDW 
5p 4 



'Pa, 



a 



1,2 



P(f) 



P3 



■) P4 



P-k/2. 



(2.24) 



Because of the sign change of 4>P3, we again, as with the 
SP phase, have a redefinition of A3 and reassociation of 
the kinks with the physical electrons. The shift of <j)pA, as 
with the S-Mott phase, only leads to phase multiplication 
of certain fields. 



E. Triality 

Of the SO(2N) groups, 5*0(8) possesses a special sym- 
metry called triality that rotates its three fundamental 
representations (the vector and the two spinors) among 
one another. It is unique to 50(8) as only in this case 
do the spinor and vector representations have the same 
dimension. The symmetry is isomorphic to Z% (hence its 
name) and so has a generator, g, such that g = 1. 

To exhibit the action of g upon the three represen- 
tations, we consider its be havio ur on the four Cartan 
bosons, 4> a p, introduced in 2.7. Under g the 



s arc 



linearly transformed as follows 



IP— * 


{<PlP + 


4>2P + 


4>ZP + 


04P)A 


2P^ 


(01P + 


4>2P - 


hp - 


04p)/2; 


3P^ 


(<PlP - 


4>2P + 


hp " 


04p)/2; 


4P^ 


(<p 1P - 


(j>2P - 


03P + 


4>ap)/2. 



(2.25) 

Given that the four Cartan bosons, <p a p are identified 
with the four quantum numbers, Ni, g acts to redefine 
them correspondingly. 

As the fermions are represented by e^ aP and the kinks 
by e ta , we can see the triality transformation acts to 
take the fermions to the kinks of even chirality (i.e. the 
number of components of a that are negative is even), 
the even kinks to the odd kinks, and the odd kinks to 
the fermions. 

We will use triality to fix some of the properties of 
the form factors in Section 5. Via triality, we can re- 
late form factors involving fermionic and kink fields and 
excitations, thus constraining them. 



III. PRESENTATION OF RESULTS 



In this section we present the computations using form 
factors of the optical conductivity, er, the single particle 
spectral function, A(k, w), and a tunneling I — V curve. 
These results are generic to any of the four massive 50(8) 
phases of the Hubbard ladders/carbon nanotubes. The 
optical conductivity is the same in each of the phases as 
they share a common 50(5) subalgebra that contains the 
conserved electric charge. The computation of the spec- 
tral function and tunneling I — V curve is done with re- 
spect to generic kinks. There is no need to specify which 
exact kink one is working with as one obtains identical 
results independent of the kink type (predominantly be- 
cause of the action of the 50(8) symmetry). Thus the 
results are independent of the identification of the kinks 
with the physical electronic excitations one makes in each 
of the massive phases. 

In the Sections 4 and 5 that follow, we lay out the 
calculation of the form factors. For those uninterested 
in the details, the bulk of these sections may be skipped 
with only the results at the end of Section 5 referenced. 



A. Optical Conductivity 



In this section we consider the response of the lad- 
der system to an electric field polarized along the legs. 
Apart from the treatment in jjj, this problem has been 
examined previously, both theoretically |34j] and experi- 
mentally j3q| . However these two latter papers did not 
consider undoped ladders at zero temperature. 

In linear response, the optical conductivity is given by 



Re[a(uj,k)] = Im[ 



A(w,k)- 



(3.1) 



where A is the current-current correlator 



A(w,fc)= / dxdTe^ r e M (T(J(x 1 T)J(0,0)))\^^ +s . 

(3.2) 



J is given by 2.20 



(3.3) 



To compute the correlator, (T(G} 2 (a;, r)G} 2 (0, 0))), we 
insert a resolution of the identity between the two J's, 
turning the correlator into a form factor sum. We then 
have 
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<T(Gi 2 (*,T)G} 2 (0,0))) = f; £ [^...^ { G?{Q)\Al{e l )...Al n {6 n )) 

n— ai,"-,a n 

n n 

x (A a „ (0 n ) • • • A ai (#i)|G} 2 (0)) exp ( - M £ m Qi cosh(^) + ^ m as sinh(^)) , 



(3.4) 



where the first sum ^ n runs over the number of par- To compute this sum in its entirety is generally an in- 

ticles in the form factor expansion and the second sum tractable problem. The usual solution is to truncate the 

Y] a . runs over the different particle types. We have also sum at some n. Here we will content ourselves with a 

extracted the spacetime dependence of each term. truncation at the two particle level: 

/df) 
^{GfmA\ 2 {9 1 )){A 12 {e 1 )\Gfm 



d6 1 d9 2 
1 2tt 1 2tt 



x exp ( — |t|\/3to cosh(6'i)) + ixV3msmh(0i)) 
exp ( — |T|m(cosh(6>i) + cosh(#2)) + ixm(sinh(9i) + sinh^))) 
(G\ 2 (0)\AU9i)Al{9 2 )){A b (6 2 )A49 1 )\Gl 2 (0^ 



a b 



(3.5) 



a/3 



The first term gives the single particle contribution to the 
correlation function. The only particle that contributes 
here is Ai 2 , denoting one of the particles belonging to 
the rank 2 tensor multiplet. At the two particle level a 
variety of contributions are non-zero. The second term 
in 3.5 gives the contribution of two Majorana fermions 



while the third term gives the contribution of kinks with 
the same chirality. 

As indicated in the introduction, this truncation of the 
form factor sum is better than it may at first seem. Be- 
cause the correlator is evaluated at zero temperature in 
a massive system, the higher order terms make contribu- 
tions only at higher energies, lu. That is, the massiveness 
of the system leads to particle thresholds. The next con- 
tribution comes from a three particle combination of even 
kink/fcrmion/odd kink that carries mass 3m. Thus for 



lu < 3m, this term gives no contribution to Re[cr(w,fc)], 
for arbitrary k. Hence our result for Re[cr(uj, k)], is exact 
for lu < 3m. 

In the case when u> does exceed 3m, we expect the 
higher particle form factors to make only a small contri- 
bution to a(u>). In practice, form factor sums have been 
found to be strongly convergent for operators in massive 
theories p^ , |2^ -|28[| . To obtain a good approximation to 
correlators involving such fields, only the first few terms 
need to be kept. Even in massless theories where there 
are no explicit thresholds, convergence is good provided 
the engineering dimension of the operator matches its 
anomalous dimension ]24|j25j ]. 

Using the results for the form factors of Section 5, we 
can put everything together and write down an expres- 
sion for Re[cr(w)]: 



-9(w 2 



Am 2 ) 



12m 2 



dx G c (x) 2 
x s(x) 



(x/3)} 



j\fbj 2 — k 2 — 4m 2 



(lu 2 — k 2 — 3m 2 ) 2 (lu 2 - k 2 fl 2 
■ i'xd / — (1 — c(x) cos(— — ))] , 

X S(X) 7T 



: exp [ / 
Jo 



(3.6) 
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where s(x) = sinh(x), c(x) — cosh(x), and 



0.5 



712 



cosh 1 [- 



k 2 - 2m 2 - 
2m 2 - 



(3.7) 



As indicated in Section 5, Aq is an arbitrary constant 
normalizing all current form- factors while G c (x) can be 
found in (5.68) . 

In Figure 1 we plot the real part of the optical conduc- 
tivity for wavevector k = 0. There is no Drude weight 
as we are at zero temperature and zero doping. We see 
there is an exciton type peak at lu — \/3m correspond- 
ing to the single particle form factor contribution. The 
first vertical dashed line marks out the beginning of the 
two particle form factor contribution to the conductiv- 
ity. The onset of the two particle contribution behaves 
as \Ju) — 2m and not as 1 /\fu) — 2m as would be expected 
in a free theory due to the divergence in the density of 
states, the van-Hove singularity, that generically occurs 
in one dimensional systems. This singularity is removed 
by the corresponding current matrix element which be- 
haves as (uj — 2m) with lu — > 2m + as is expected generi- 
cally because the low energy behaviour becomes strongly 
renormalized in the presence of even weak interactions 
f§- 

The optical conductivity was computed in (lj using the 
large N limit of SO(2N), or in an alternate language, 
an RPA approximation. In such an approximation, the 
model becomes equivalent to a theory of four massive, 
non-interacting Dirac fermions. Hence finds that the 
van-Hove singularity is present. 

The second vertical dashed line in Figure 1 at 3m 
marks the point where the three particle form factors 
would begin to make a contribution. Up to this point, 
the result is exact. We note that the three particle contri- 
bution is strictly a consequence of interactions. In 5*0(8) 
language, two kinks of opposite chirality together with a 
fermion will couple to the current operator. In a free the- 
ory these different particles would not all exist and there 
would be no three particle contribution. 

If we were to compute the three particle contribution, 
three possibilities present themselves. The three particle 
density of states approaches a constant as lu — > 3m + . If 
the corresponding matrix clement vanishes as lu — >■ 3m + , 
the contribution opens up gradually, leaving a (to) contin- 
uous at to = 3m. If the three particle matrix element also 
approaches a constant value as to — > 3m + , the conduc- 
tivity is marked by a jump at uj = 3m. But if the matrix 
element diverges in this limit, we expect a corresponding 
divergence in the conductivity at to = 3m. 



0.4 



0.3 



0.2 




Figure 1: Plot of the optical conductivity at 
wavevector k = 0. 



In Figure 2 we plot the real part of the optical conduc- 
tivity for wavevector k = ir/3. As a result of the form of 
Re[(cr)] (see 3.6), shifting k from to ir/3 shifts the single 
particle form factor contribution to m\J tt 2 /9 + 3 and the 



onset of the two particle contribution to my/n 2 /9 + 4. 
We see that as k grows, the contribution from the single 
particle moves closer to the two particle threshold. 

It can now be asked how perturbations to SO (8) Gross- 
Neveu will affect the computation of the optical conduc- 
tivity. We consider this in the broadest terms by focusing 
upon how the spectrum of 50(8) Gross-Neveu is changed 
under a perturbing term. We do so through straightfor- 
ward stationary perturbation theory, in the same spirit 
that [|TJ treated the off-critical Ising model in a magnetic 
field. The most general possible perturbation takes the 
form 



-ffpert — XG a bG 



(3.8) 



where G a b, G c d are SO (8) currents (of unspecified chiral- 
ity). For such a perturbation it is necessary to consider 
degenerate perturbation theory. Thus in a given parti- 
cle multiplet (for example, the fundamental fermions in 
the vector representation), the perturbed energies arise 
through diagonalizing the matrix 



Ma = (A]{d)H v 



(3.9) 



where here the index i,j indicates the particles A^Aj 
belong to the multiplet of concern. In the case that 
G a b = G c d, Mij is necessarily diagonal, i.e. nondegener- 
ate perturbation theory is sufficient. 
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fashion. Like the unperturbed case, there will be sev- 
eral two particle contributions to the optical conductiv- 
ity. However unlike the unperturbed case, the thresh- 
olds of the two particle contributions will not all oc- 
cur at u> = 2m but be distributed about this energy. 
Thus the two particle contribution is arrived at (approx- 
imately) by superimposing several slightly shifted two- 
particle contributions similar to those found in Figures 
1 and 2. But given the optical conductivity vanishes at 
threshold, the qualitative picture remains effectively un- 
changed (i.e. the superimposed contributions will appear 
nearly identical to the original picture). That the opti- 
cal conductivity vanishes at the two particle threshold is 
a result of the vanishing of the relevant matrix element 
at threshold. This should be robust under perturbation 
as it is ultimately a consequence of the mere presence of 
interactions and not some particular type of interactions. 



Figure 2: Plot of the optical conductivity at 
wavevector k = tv/3. This plot is exact up to 
w/m = -\/ir 2 /9 + 9 at which point the three par- 
ticle states begin to make a contribution. 

It is important to emphasize that this procedure can 
be han dled in the context of integrability. The expression 
in 3.9 is no more than a form- factor which can readily be 
computed. Moreover as the theory is massive, pertur- 
bation theory is well controlled. We expect the unper- 
turbed theory to describe all qualitative features of the 
model while the perturbations to only introduce small 
quantitative changes. 

The consequences of a perturbation for the optical 
conductivity are two-fold. We expect the exciton peak 
(found, for example, in Figure 1 at k = and ui = \/3) 
to split. In the unperturbed model the peak results from 
a single rank two bosonic bound state coupling to the 
current operator. When the matrix M a b above is diago- 
nalized, this particular state should be mixed into many 
others resulting in several states that couple to the cur- 
rent operator. However we do not expect the functional 
forms of the exciton peaks to change: they should remain 
delta functions. They must do so provided the pertur- 
bation is not so large as to push the exciton peak past 
the threshold of two particle states where it then con- 
ceivably could decay. As there is a gap between the ex- 
citonic peak and the two particle threshold, this will not 
happen for small perturbations. Experimentally changes 
to the exciton peak may not be detectable. Given that 
any experiment will be conducted at finite temperature, 
the excitonic peak will be thermally broadened, perhaps 
washing out any splitting of the original zero temperature 
peak. 

We also expect the perturbation to affect the onset of 
the two particle threshold, although in a less dramatic 



B. Single Particle Spectral Function 

In this section we compute the single-particle spectral 
function of the electrons of the Hubbard ladder. To do 
so we first consider the correlator, 



G(k x ,k y ,T) = 



1 = 1,2 



dx e 

x<T(o, o (a;,T)o, t a (0,0))). (3.10) 



Here k y takes on the values 0, 7T. We then define the 
particle/hole spectral functions, A p / hl as follows: 



A p (k,w) +A h (-k,-Lj) 
Im 



dr e'^ T G(k x ,k y ,T)\ 



. (3.11) 



In keeping with |5J], we have not explicitly summed over 
spin, a. 

As described in Section 2, electronic excitations around 
the Fermi point correspond in the Gross-Neveu language 
to low energy excitations of kinks. We thus expect to 
recast the Greens function, G, above in terms of kink 
correlators. This in fact can be done with the result, 

G(Pk Fi + k, k y u t) = 

dxe ikx (T(c Pla (x, r)4fa (0, 0))) , (3.12) 

where i = 1,2 and k y i = (2 — iW. The c's, the bonding- 



anti-bonding electron s give n in 2/2, are in turn related to 
the various kinks via 2.18. The Greens function on the 



r.h.s. of 3.12 is thus equal to 

<T(c Pm (x,r)4 M (0,0))) = (T( Ka V>±(*,TKV>S(0,0))>. 

(3.13) 
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where a (a being its charge conjugate) is the partic- 
ular kink corresponding to the Fermi point {kpi,k y i). 
The K a are Klein factors included to ensure the ip a 
are anti-commuting. Because of the 50(8) sym- 
metry together with its associated triality symmetry, 
(T(/c Q ,V±( a; , t)k s ^^(0, 0))) turns out to be independent 
of the type a of kink. It is only sensitive to whether the 
kink field is right (+) or left (— ) moving. 

To compute this correlator, we again expand to the 
two lowest contributions: 

(T(/s a ^(a: ) r>0)/s a ^(0,0))) = 

^W(^r)4(<?i))(^«(e l )V?(0)) 
+ E f y^W(*.r)4(*)4(*)> 

afj ^~ oc 77 7F 



x (4,(01)4,(02)^(0)). (3.14) 



The first contribution, the one particle contribution, 
comes from the kink excitation, A & , destroyed by the 
field, ift a . The second contribution, a two particle contri- 
bution, arises from kinks, Ap, of opposite chirality to A a , 
and Majorana fermions, A a . (This reflects the group the- 
oretical fact that the tensor product of an SO (8) spinor 
representation with an 5*0(8) vector representation gives 
the other SO (8) spinor representation j3f|.) The first 
contribution not included is a bound state-kink pair. It 
begins to contribute at u = (1 + v3)m. 

From the form factor expressions from Section 5, we 
can then write down the expression for the spectral func- 
tions, A p / h (uj,k), 



A p (uj, Pkpi + k, kyi) = Ah{u>, -Pkpi + k, ky%) = 



7T \Cp\ 



Pk 



+9(uj - \Jk 2 + 4m 2 ) 
x exp [ 



TO 

16m 4 Al 



Vk 2 



l 



\6(w — \J k 2 + m 2 ) 
1 



UJ-Pk (w 2 

dx Gf(x 



o x s(x) 



2 _ m 2~)2 ^Jl 

(1 - c(x) cos(-^.))], 



k 2 — 4m 2 



(3.15) 



where Ap is the (unspecified) normalization of the two 
particle kink form-factor, Gf(x) is given in (5.70), and 
c± is found in (5.73). For P = R = + (i.e. right-moving 
electrons/kinks), this function is plotted in Figure 3. 




Figure 3: Plot of the single particle spectral func- 
tion for right moving kinks. The more darkly 
shaded region corresponds to greater spectral 
weight. 



The parabolic line in Figure 3 arises from the sin- 
gle particle form factor contribution, and represents the 
standard dispersion relation of a particle of mass, to. 
Above this curve comes the two particle form factor con- 
tribution to the spectral function. This contribution is 
bounded by the curve, u> = \/k 2 + 4m 2 , and so the sin- 
gle particle states do not cros s into the two particle re- 
gion. As can be seen from 3.15 , the two particle contribu- 
tion opens up at threshold with a square-root singularity, 
indicative of the van-Hove singularity in the density of 
states. 

The plot is manifestly chiral with weighting greater 
for k > than for k < 0. This is to be expected as we 
are plotting the excitations linearized about the Fermi 
momentum, +kpi. The heavier weighting for k > indi- 
cates that is easier to create excitations above the Fermi 
sea than below it. It is interesting indeed that excita- 
tions below the Fermi surface can be created at all and 
is a mark that interactions are at play. 



As 3.15 indicates, we have computed the single par- 
ticle spectral function in the neighborhood of the Fermi 
points. However in JjJ], the authors found two particle 
spectral weight to exist away from the Fermi points at 
harmonics k — ±3kpi and uj = 2m, that is the same 
energy threshold as at the Fermi points, kpi, them- 
selves. They identified the spectral weight at these points 
through combining the SO (8) picture of the ladders with 
the original lattice formulation. As such, we are unable 
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to say anything definitive about spectral weight at these 
points as we are forced to remain with the 50(8) descrip- 
tion. 

The effect of a perturbation breaking 50(8) can be 
sketched in a similar fashion to that of the optical con- 
ductivity. In the case of the single particle contribution 
to the spectral function, we expect that under degenerate 
perturbation theory the single parabolic line in Figure 3 
to split akin to the behaviour of the excitonic peak in the 
optical conductivity. Again this effect will be masked by 
finite temperatures which thermally broaden such fea- 
tures. 

At the two particle level, both in the unperturbed and 
perturbed cases, several sets of two-particle pairs con- 
tribute. In the unperturbed case, these sets all begin to 
contribute at oj = \J Am 2 + k 2 . In the perturbed case, 
the sets will begin to contribute at different points. For 
small shifts, a plot such as that in Figure 3 will remain 
much the same. However it is worthwhile to point out 
that with each two-particle pair is associated a square- 
root divergence due to the van-Hove singularity in the 
density of states. Thus we would expect to find a series 
of such singularities in the presence of a perturbation. 



C. Tunneling Current 

In this section we study the tunneling between a 
metallic lead and the carbon nanotube/Hubbard ladder 
through a point contact. Our starting point is a La- 
grangian describing the nanotube/ladder, the metallic 
lead, and the tunneling interaction: 

£ = AsO(8) + £lead + Ctun- (3.16) 

£-SO(8) is the Lagrangian of the 50(8) Gross-Neveu 
model. 

The electron gas in the lead is, in general, three di- 
mensional. However, in the context of tunneling through 
a point contact, the electron gas can be mapped onto an 
one dimensional chiral fermion (see for example fl37j-]39| ) . 
The general idea is well illustrated by its application to 
the Kondo problem. There an electron scatters off a spin 
impurity at x = 0. The scattering is determined by the 
electron operator, ip(x = 0). As ^(0) only depends on 
its spherically symmetric, L = 0, mode, one can consider 
the scattering electron in terms of an ingoing and out- 
going radial model defined on the half-line, r £ [0, oo]. 
Unfolding the system onto the full line leaves one with 
a chiral fermion. Wc emphasize however that the map 
requires no special symmetry; the result is exact regard- 
less of particular anisotropies (39). As a consequence, we 
write Ci eadl as 

Ulead = (3.17) 



where is a massless, left moving fermion, and z = 
(r + ix)/2. 

It remains to specify Ctun- In order to preserve charge, 
the electrons must couple to the kinks of the 50(8) 
Gross-Neveu model, the excitations with the quantum 
numbers of the electron. Thus 

+flfl[* t (r)^(r)+V?(r)*(r)]<y(a:). (3.18) 

Here we have coupled the lead electrons to both the right 
and left moving fields creating the kink, a, and have al- 
lowed the two couplings, and gn, to be unequal. How- 
ever as we will work to lowest non- vanishing order in the 
tunneling matrix elements <?_l/_r, the tunneling current 
will depend upon the sum, g\ + g^, that is, the contri- 
bution of the left and right channels to tunneling will 
add linearly Similarly, permitting other kinks to couple 
to the lead electrons will give lowest order contributions 
which simply add. 
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Figure 4: Plot of the tunneling current as a func- 
tion of applied voltage. The dashed curve de- 
scribes tunneling into a non-interacting fcrmionic 
system of mass, m. The vertical dashed line marks 
where a second set of two particle states begins to 
make a contribution. 

The calculation of the current to lowest non- vanishing 
order in g^m follows the standard route. The tunneling 
current operator is given by 

I(t) = is«(^(r)*(r) - 3%)V?(r)) 

+»Si(^?(r)*(T) - * f (T)^(r)). (3.19) 
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In order to induce current flow, (I(t)), through the point 
contact, one biases the lead with a voltage, V. This bias 
can be taken into account via a gauge transformation, 



*(r) -» e lVT #(r) 



(3.20) 



In effect we have shifted the energy levels of the electrons. 
Treating the couplings, gh/gn, with linear response the- 
ory, we find 

(I(uj)) = g 2 L Re\ f drie^l 



e- iyT (C(r)C(0))(*t( T )*(0)) 
-e WT (C (r)C (0)) (*(t)*+ (0))] | v -^+ 5 



(3.21) 

The lead electron correlator (^(t^^)) is well known: 

1 



(&(t,x)V(0)) = (*(t, 1)^(0)) = 



(3.22) 



With this it is straightforward to express the dc current, 
(J(cj = 0)), in terms of the single particle kink spectral 
function, 

(J(w = 0)> = — / doj dk[glA^w,k) 

+g 2 R A + (uj,k)], (3.23) 

where A±(w, fc) = A p (u>, ±kpi + k, k y i) + Ah{— ui, T ^Fi — 
k, kyi), and A p / h are the spectral functions given in [3.15 



We note that as a technical point, in deriving the above 
equation we have displaced, 'J', the lead electron opera- 
tor, slightly from x = 0. In this way we cure the UV 
divergence attendant as r — > 0. At the end of the calcu- 
lation we then take x to 0. 

In the previous section we have computed A±(to,k) 
exactly for energies u) < (\/3 + l)m. Inserting |2.15| into 
3.23, we find (1(0)) takes the form, 



(HO)) 



— (fffl + AW 2 - m 2 ) l / 2 9(\V\ - m)sgn(V) 



+0(|V| — 2m)sgn(V) x two particle contribution. 

(3.24) 

We see that for | V\ < 2m, the system behaves as a gapped 
free fermion. The first sign that there is any interaction 
comes for \V\ > 2m where the voltage begins to probe the 
two particle states, a signature of interacting fermions. 

We explicitly plot (1(0)) in Figure 4. The square root 
behavior near V/m = 1 and subsequent linear form is 



typical of a gapped fermion. At V/m = 2, the two parti- 
cle states begin to contribute leading to a small change in 
the slope of the I —V curve. At V/ m = v3 + 1 (marked 
by the the vertical dashed line), a second set of two parti- 
cle states (a mass \/3m bound state together with a kink) 
begin to contribute and at this point the result ceases to 
be exact. However as with the current correlators, we 
expect this higher energy contribution to be small. 

The change in slope in the I — V curve at V/m — 2 can 
be explicitly computed. To do so we consider dv(I(0)). 
This quantity is given by 



d v (I(0)) = -(«? 

7T 



g 2 L ) dkA ± (V,k). 



(3.25) 



We can thus see dy(I(0)) directly measures the local den- 
sity of states at x = of the nanotube/ladder system. 

We plot dy(I(0)) in Figure 5. The square root singu- 
larity at V = m signals the singularity of the density of 
states in an one dimensional system. At V = 2m we see 
a sudden jump, indicative of the onset of the two particle 
contribution. The height of the jump can be determined 
exactly: 

d v (I(0))(V/m = 2+) - d v (I(0))(V/m = 2~) = 
8A 2 F/ 2 2N r f°° dxG f (x) , , ... . 

The region m < V < 2m of dy(I(0)) completely deter- 
mines m (by the location of the jump), as well as an over- 
all scale (the product of (g\ + g\) and the constant Ap, 
normalizing the spectral function). Dividing out these 
non-universal numbers leaves a universal number, char- 
acterizing the magnitude of the jump: 



8 .r 

9 exp[ y 



dx Gf (x) 
x s(x) 



(l-c(x))}. 



(3.27) 



This number represents a definite prediction based upon 
the integrability of the model. 

We now consider the approximate effect of pertur- 
bations breaking integrability on the tunneling conduc- 
tance. As the tunneling conductance is determined di- 
rectly from the single particle spectral function, we can 
deduce how the former is affected from how the latter 
is changed. At a given momentum, the single particle 
contribution to the single particle spectral function un- 
der perturbation comes at a discrete set of energies. In 
terms of the tunneling conductance, we expect a series 
of closely spaced square root divergences (a sawtooth be- 
haviour) about V — m indicative of a series of van-Hove 
singularities. As the perturbation is removed these sin- 
gularities would collapse on top of one another leaving 
the original picture in Figure 5. 

In the unperturbed case the two-particle threshold is 
characterized by a jump in the differential conductance. 



1G 



Under a perturbation, this jump would become a stair- 
case or a series of smaller, closely spaced jumps. This is 
a reflection of the series of van-Hove singularities found 
about lo = \/k 2 + 4m 2 in the two particle contribution 
to the single particle spectral function. 
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Figure 5: Plot of the differential conductance as 
a function of applied voltage. The dashed curve 
marks the single particle contribution to this quan- 
tity while the solid curve gives both the single and 
first two particle contribution. The latter plot 
is exact up to V/ m — V3 + 1 (indicated by the 
dashed vertical line) where a bound state-kink pair 
begins to make a contribution. 

IV. S-MATRICES FOR 50(8) GROSS-NEVEU 

In this section we review the S-matrices for scattering 
among the kinks and elementary (Majorana) fermions. 
From these we derive S-matrices for scattering between 
Dirac fermions, kinks and Dirac fermions, and rank 2 
bound states and Dirac fermions. Knowledge of these 
S-matrices will allow us to compute the various form fac- 
tors needed to obtain the results in the previous section. 
Before beginning it is necessary to lay out the group the- 
ory of 5*0(8). In particular it is necessary to describe the 
spinor representations (i.e. the kinks) of 50(8). 

A. 5*0(8) Spinor Representations 

The spinor representations are expressed in terms of 
50(8) 7-matrices. The 7-matrices in turn are built out 
of the two-dimensional Pauli matrices, a l, s. As we are 
interested in 50(8), we consider four copies of the er"s, 



a\,al,ol,<y\, (4.1) 

each acting on a different two-dimensional space: 

a l n \a>i,--- ,a 4 ) = \a x , ■ ■ ■ , a l n a n , ■ ■ ■ , a 4 ). (4.2) 

A particular basis vector, |a), in the corresponding 16 di- 
mensional vector space can be labeled by a series of four 
±1/2, i.e. 

a = I ±1/2, ±1/2, ±1/2, ±1/2). (4.3) 

Physically, the kink associated with \a) then carries ±1/2 
of the four U(l) quantum numbers, Ni, corresponding to 
the Cartan elements of 50(8), discussed above. 

In terms of the ct"s, the 7-matrices are defined by (fol- 
lowing the conventions of |40|l), 

n-l 

72n-l Y[vh 
n-1 

J2n=<jl®n<jl; (4.4) 

for 1 < n < 4. These matrices satisfy the necessary 
Clifford algebra: 

{ 7a , 7 fc} = 2(U. (4.5) 

In this representation the Clifford-algebra generators, 7„, 
are imaginary and anti-symmetric for n even, while for n 
odd, they are real and symmetric. The 50(8)-generators 
are represented by 

a ab = 1( 7 V-7V), (4.6) 

in analogy with the more familiar 50(4) case. 

The 16-dimensional space of the 7's decomposes into 
two 8-dimensional spaces, each of which forms one of the 
two irreducible 50(8) spinor representations. The de- 
composition is achieved explicitly via the hermitian chi- 
rality operator, T, 

8 

r = Y[j a . (4.7) 

a=l 

F is such that it commutes with all 50(8) generators and 
is diagonal with eigenvalues ±1. If la^) is a state with an 
even (odd) number of negative components (i.e. states 
with positive or negative isotopic chirality), then 

r|a ± )=±|a ± ). (4.8) 

Thus the operators, (1 ± F)/2, project onto the two irre- 
ducible subspaces. 



17 



The last item to be presented in this section is the 
charge conjugation matrix, C. In terms of the 7's, C is 
given by 



C = 72747678- 



(4.9) 



C is completely off-diagonal (as expected). If \a) is the 
anti-particle of a kink, \a), then 

\a)=C af3 \/3). (4.10) 

C is such that C 2 = 1, C T = C, and CYC" 1 = 7 a *. 



B. S-Matrices for the Elementary Fermions and 
Kinks 



To describe scattering we introduce Faddeev- 
Zamolodchikov operators, A a {8) and A a (8), that cre- 
ate the elementary fermions and kinks respectively. 8 is 
the rapidity which parameterizes a particle's energy and 
momentum: 



p = msinh(#); E = mcosh(9). 



(4.11) 



Because the SO (8) Gross-Neveu model is integrable, 
scattering is completely encoded in the two-body S- 
matrix. This S-matrix, in turn, is encoded in the commu- 
tation relations of the Faddeev-Zamolodchikov operators: 

Ai{e 1 )A 2 {9 2 ) = Sf%{8 1 2)A 4 (6 2 )A 3 (6 1 ). (4.12) 

^12(^12) is the amplitude of a process by which particles 
1 and 2 (be they kinks, fundamental fermions, or bosonic 
bound states) scatter into 3 and 4. It is a function of 
812 = 81—82 by reason of Lorentz invariance. 

Scattering between the fundamental Majorana 
fermions, first determined in EL is given by 

S c a d b (0) = & ab 5cdOi{8) + S ac 6 bd a 2 (8) + S ad S bc a 3 {8), (4.13) 

where the cr's are defined by 

6 ITT — 8 

^(0) = -i^a 2 (6); 

a 2 {9) = -Q{d)Q{m-9) s{ fl +i ^f (4.14) 
s(8) — iv3/2 

where s{8) = sinh(6>). Q is given to be 

nm r(l/6-ig/27r)r(l/2-jg/27r) 

Q{6) = r{-iO/2*)r(2/3-M/2*) " (4 ' 15) 

S%h(9) satisfies a Yang-Baxter equation, 

S££ (^12)^(^3)^(^23) = 



SZZ {^szli (6 13 )S b c l b c > (8 12 ). 



(4.16) 



Physically, the Yang-Baxter equation encodes the equiv- 
alence of different ways of representing three-body in- 
teractions in terms of two-body amplitudes. Formally, it 
expresses the associativity of the Faddeev-Zamolodchikov 
algebra. The S-matrix, S%£(6), also satisfies both a cross- 
ing, 



s%(e) = s£(i*-o) 



and unitarity relation, 



S c a d b (0)S e J(-8) = S e J f b . 



(4.17) 



(4.18) 



These constraints determine S^f up to a 'CDD'-factor. 
Such factors determine the pole structure in the physical 
strip, Re(0) = 0, < Im6> < 2ir, of the scattering matrix 
and so are indicative of bound states. Here the CDD 
factor is 



s(8) + iV3/2 
s(0) - n/3/2' 



(4.19) 



and reflects the fact that two fermions can form a bound 
state of mass ^/3m transforming as a component of the 
29 dimensional representation of the 5*0(8) Yangian. As 
stated in the introduction, this representation is formed 
out of a rank-2 anti-symmetric 5*0(8) tensor and a scalar. 
The overall sign of the S-matrix is determined by exam- 
ining the residue of the pole in 5°° at 8 = iir/3. This 
pole is indicative of the formation of a mass y/im scalar 
bound state in the s-channel and so should have positive 
imaginary residue. 

Scattering between the kinks and fundamental 
fermions is given by 



Sli = -h{e)8 ab 8 otf) -t 2 {e)af a , 



The form of 4.18. first written down in 



(4.20) 



is determined 



by writing down all possible 50(8) invariant tensors with 
indices a, b and a, j3. t\ and t 2 are described by 



h(8)-t 2 (8) 
h(8) + 7t 2 (8) 



8-i2ir/3 
8 + i2n/3 



(4.21) 



where 82(8) is given by 

s{8) +i v / 3/2r(5/6- 



5*2(0) 



■i6/2ir)r(4/S-i6/2n) 



8(0) - iV3/2 T(5/6 - »0/27r)r(4/3 - 



■i0/2n) 
(4.22) 



Again, the Yang-Baxter equation, crossing, and unitarity 

These constraints 



are used to determine 4.21 



and 4.22 



were first explicitly solved in [J40[ . The overall sign is 
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fixed by insisting the pole in 5"^ at 9 = i2ir/3, indica- 
tive of a kink formed in the s-channel as a fermion-kink 
bound state, has positive imaginary residue. 
Kink-kink scattering takes the form 



Q -y5 _ J_ V"^ U r {9) (r) (r) 
r— 



(4.23) 



where is a rank-r anti-symmetric tensor: 

Here A represents a complete antisymmetrization of the 

(r) (r) 

gamma matrices. By cr^a s J we mean a trace over all 
possible rank-r anti-symmetric tensors: 

(r) (r) \ ai---a r ai---a T (a 

Again the generic form of was determined in [42|, 
while the specific forms of the it's were given in [40|. 
There it was found 

u i+r {9) = (-iy Ui - r (9); 

(tt\ , a , (l-r/3)-(l + z9/ir) . 

Ur+2{9) = ^(i-r/aj + a + W 

M 2 (^) = CT2(«7T - 0) - <7 3 (iTT - 9)\ 

ux{0) = S 2 {m - 6). (4.26) 

The 5*0(8) Gross-Neveu model has an isotopic chiral- 
ity conservation law ffl]. Thus opposite chirality kink 
scattering is determined solely by u n , n odd, while same 
chirality kink scattering is determined solely by u n , n 
even. The overall sign of the S-matrix again is chosen 
so that appropriate s-channel poles have positive imag- 
inary residues. That the it's bear a close relationship 
to amplitudes from fermion-fermion scattering and kink- 
fcrmion scattering is a reflection of the triality symmetry 
of SO(8). 



C. S-Matrices for Dirac Fermions 

On occasion the Majorana fermion basis will not be the 
most convenient. For example, as these fermions do not 
carry definite U(l) charge, they will not couple nicely to a 
chemical potential introduced when the system is doped. 
Thus we would want to consider, as in |^|, a fermion ba- 
sis with well-defined U(l) charges, i.e. we would want to 
reexpress the Majorana fermions as Dirac-fermions. 

As there are eight Majorana fermions in 5*0(8) Gross- 
Neveu, we have four Dirac fermions. Denoting the cor- 
responding Faddeev-Zamolodchikov operators by Df(9), 
where D~(9) is £>+(#) 's anti-particle, we have 



Df{6) = -j={A 2i {6)±iA2i-x{6)) 1< ? <4, (4.27) 

where the Ai's are Majorana fermions. Because scatter- 
ing is determined by the commutation relations, 



h)D^{9 2 ) = S\ 



(fec fc )(iC() jjc k 



iCi){jCj) ~ k w)D?Wi), (4-28) 



where c = ±1 is indicative of the U(l) charge carried, we 
can 
ing, 



can determine Si ko \}. 1 ) from |4.13|. We find the follow- 



+SuS CzCl S :jk S CkC:j a 3 {9), 



(4.29) 



where Cp = o\ is the charge conjugation matrix for Dirac 
fermions. We point out that two Dirac fermions of the 
same U(l) charge will scatter diagonally, i.e. 

Df{9 l )Df{9 2 ) = (a 2 (9 12 ) + a 3 (9 12 ))Dt(9 2 )Dt(9 1 ). 

(4.30) 

We will also want to con sider scattering between kinks 
and Dirac fermions. From 4.20 we find 



^{Yir_ Ci ) 0a t 2 (9). (4.31) 



where 



: 7 



23-1 ± rfl 



(4.32) 



When A a {9) carries U(l) charge Cj/2 (i.e. when the U(l) 
charges of A a {9) and Di carry the same sign) the above 
form simplifies greatly: 



(4.33) 



i.e. scattering becomes diagonal. This is exploited in M 
to compute the excitation energy of a kink in the doped 
system. 

We want to consider one last S-matrix involving the 
Dirac fermions, that of a Dirac fermion with rank 2 tenso- 
rial bound states. We will not do this in general. Rather 
we consider a bound state carrying two non-zero positive 
charges, say Ui = +1 and Uj = +1, scattering off a Dirac 
fermion carrying charge either Ui = +1 or Uj = +1. This 
is precisely the situation encountered when treating spin 
1 excitations in the doped system. 

We thus introduce the bound state Faddeev- 
Zamolodchikov operator, As (9). As (9) can be thought 
of as a bound state of two Dirac fermions with charge 
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Ui = +1 and Uj — +1, and consequently can be rep- 
resented in terms of the corresponding two Faddeev- 
Zamolodchikov operators, 

i 9(i+)(j+) A B(d) = 

res^o A i+ (6 + S + iu^ s )A j+ (9 - iufV^), (4-34) 



U+)B> 



where resa = o denotes the residue at 8 — and B is the 
charge conjugate of B. In writing 4.34 , we have taken the 
particle normalization to be (6\8') = 2ir5(8 — 6'). The u's 
are given by 



u sB = * - u 



sBi 



(4.35) 



where iu s sh is the location of the pole indicative of the 
particle s 1 in the s—B S-matrix. Here we have 



\i+)B 



a 



(i-) 



5tt/6. 



(4.36) 



The g's are related to the residues of the poles in the 
(£+)— (i+) S-matrix and can be interpreted as the am- 
plitude to form the bound state from (i+) and 

fi+)U+) * s ^ e l° ca ^i° n °f the P°-^ e indicative of B, 
is defined by 



If u B 



s, 



(0) 



Q B a (i+)ti+) 
■ g (»+)(j+) g g 

6 ~ lU (i+)U+) 



(4.37) 



It is then easy to show (from 4.30) that 



9{i+){j+)9B 

n B Ji+)(J+) 



So', 
—So; 



(4.38) 



where So is some constant. Hence gE + \r i+ \J 'gf i+ \tj + \ = 
— I. This last relation is all we will need to determine the 
scattering between Ag(9) and Df-. 



From |4.32| and |D| we find that SgW,^ is given by 



S 



B(i/i+) 



B(i/i+)W = a2( - 9 ~ i7r / 6 )(°a( e + in/6)+a 3 {8 + in/6)) 
-a 3 (9-m/6)a 3 {9 + i%/6) 



x (a 2 (9 - in/6) + a 3 (8 - in/6)), 



(4.39) 



so that scattering between Ab(0) and Df. is diagonal 



V. FORM FACTORS 

Here we determine the needed form factors to com- 
pute correlators in the undoped system. If the reader 
is uninterested in the actual derivation, the results are 
summarized at the end of the section. 



A. Basic Properties: Two Particle Form Factors 

The two particle form factors of a field ip(x) are defined 
as the matrix elements, 



ft 2 {9 1 ,9 2 ) = {i>{0)A 2 {9 2 )A 1 {9 1 )). 



(5.1) 



The form of fx 2 (9i,9 2 ) is constrained by integrability, 
braiding relations, Lorentz invariancc, and hcrmiticity. 

The constraint coming from integrability arises from 
the scattering of Faddeev-Zamolodchikov operators. For 



5.1 to be consistent with 4.12, we must have 



fUOM = $klfin)f&[fiM- (5.2) 



The second constraint can be thought of as a periodicity 
axiom. It reads 



f^(9 2 ,9 1 ) = R^ff 2 (8 1 -2m,9 2 ) 



(5.3) 



R^> 2 tI> is a phase 8 that arises from the non-locality of the 
fields in SO (8) Gross- Neveu. This non-locality leads to 
non-trivial braiding relations between fields which R^ 2 ^ 
encodes: 

ip2(x,t)ip(y,t) = R^^{y,t)^ 2 {x,t); x<y. (5.4) 

ip2 here is the field that is associated with the particle, 
A 2 (9 2 ). x_2 can be derived through crossing symmetry 
constraints (see, for example, |3^p2]] ). 

The form factor, ff 2 , must also satisfy Lorentz invari- 
ance. If tp has Lorentz spin, s, ff 2 will take the form (at 
least in the cases at hand), 



f&(0 1 ,6 a ) = e 



s(e 1 +e 2 )/2 



(5.5) 



where by virtue of #12 = 8\ — 8 2 , / is a Lorentz scalar. 



The constraints |5.2| , |5.3| , and |5.5| , do not uniquely spec- 
ify ff 2 . It is easily seen that if ~ff 2 satisfies these axioms 
then so does /^i?(cosh(0i2)), where R{x) is some ratio- 
nal expression. The strategy is then as follows. One first 



In certain circumstances J? is a matrix, thus marking out 
the difference between abelian and non-abelian braiding. 
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determines the minimal solution to the constraints, mini- 
mal in the sense that is has the minimum number of zeros 
and poles in the physical strip, Re(0) = 0, < Im(8) < 
2tt. Then one adds poles according to the theory's bound 
state structure. If the S-matrix element scattering parti- 
cles 1,2 has a pole at 6 = iu, then ff 2 (9i, 6 2 ) has a pole 
at 



y 2 = tfi + iu. 



(5.6) 



Insisting that ff 2 has such poles and only such poles fixes 
R(x) up to a constant. 

The phase of this constant can be readily determined. 
Appealing to hermiticity gives us 

(V(o)A 2 (0 2 Mi(0i)>* = (4(004(02)^(0)) 

= (V f (Q)A l {9 1 - in) A 2 (6 2 -in)), (5.7) 
where the last line follows by crossing and so 



ff 2 (Oi,e 2 r = fl(0 2 -in,6i-iir). 



(5.8) 



For the form factors we will examine, this will be enough 
to fix their overall phase. 



B. Basic Properties: One Particle Form Factors 

One particle form factors are in a sense trivial; Lorentz 
invariance completely determines their form. If ip(x) has 
Lorentz spin, s, then 



ff(6) = mo)M0)) = 



(5.9) 



where c is some constant. To determine c we use the 
theory's bound state structure. 
In analogy with 4.34 we write 



ig^A^O) = res,5 = o A 2 (6 + S + iu 3 2l )A 3 (e - iu 2 3l ). (5.10) 

Then we have 

iglsffW = resfco f&(0 - «&, 9 + 5 + ififj). (5.11) 

Thus knowledge of the two particle form factors com- 
pletely determines their one-particle counterparts. 

C. Braiding of the Fields 



In order to employ the periodicity axiom 5.2, we need 
to specify the braiding of the fields. In order to do this for 
the kink and fermion fields, we identify these fields with 
their corresponding excitations, A a and A a . The braid- 
ing of the fields is then encoded in the asymptotic lim- 
its of the corresponding S-matrices (see [p2[). Precisely, 



if pi?i2 is defined by (where P = ± denotes right/left 
movers) 

tpip(x, t)ip 2P (y, t) = P R 12 ip 2P (y, t)ip 1P (x, t), x < y, 

(5.12) 

then pRi 2 is given by 

P R l2 = (5 1 1 |(+oo)) p . (5.13) 

In this way we find 

P _ J — 1, a,/3 are the same chirality ,g 
± a/3 I ±2 5 a,/3 are of opposite chirality; 

and 

±R aP = ±i. (5.15) 

The braiding relationship among the kinks differ from 
what one finds by defining the kink fields as vertex oper- 
ators, 

V'oP = e M t (5.16) 
and using the easily derived braiding relations: 

e ia<t> P (x,t) e ib4> P {y,t) _ & iPTrab & ib(j)p{y,t) e ia<f>p(x,t) X > y 

(5.17) 

Similarly, the correct braiding of kinks with the Majorana 
fermions is not obtained by considering the braiding of 
the corresponding vertex operators. 

The only other fields that we concern ourselves with 
in this paper are the currents. However the current are 
local fields and so have trivial braiding with the kinks 
and the fermions. 



D. Two Particle Form Factors for the 5*0(8) Currents 

There are two possible two particle form factors for the 
50(8) currents, one involving two fermions, 

„/#(0i,0 2 ) = (G*(0)Mh)M6i)h (5-18) 
and one involving two kinks, 

„/$(*!, 6 2 ) = {Gf{Q)Ap{d 2 )A a {9 1 )), (5.19) 

where the two kinks have the same chirality. By the trial- 
ity symmetry of SO(8), these two form factors will have 
the same functional form. Indeed we will use triality to 
set their relative normalization, crucial to the calculation 
of the conductivity in Section 3. 

Calculation of 

To determine ^f"^, one first fixes its group theoreti- 
cal structure. Given that G^ b is anti-symmetric and that 
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A d A c must be anti-symmetric in d and c in order to cou- 
ple to G° fc |3(| , we must look for an invariant tensor 
anti-symmetric in the pairs a, b and c, d. The only choice 
is 5 ac Sbd - SadSbc, and so takes the form, 



p.fcd (01 J 02 ) — (^ac^bd — 5 a d5bc)f^ (01 , 02 ) 



Lorentz invariance demands 



/m(0i,0 2 ) 



_ (J6!+e 3 )/2 



i-ire 



-(9i+e 2 )/2 



(5.20) 

)/(012). 

(5.21) 



Hav ing so c onstrained the form of ^/"j 02), conditions 
5.2 and 5J3 tell us f(8i2) must satisfy 



/(-0) = -/(0)M0)-<7 a (0)); 
f(-0) = f(0-2ni). 



(5.22) 



As the current is bosonic, the braiding here is trivial. 
Because <J2 — 03 can be written in the form, 



s(—)G c (x) 



X ITT 



o 2 (0) ^ C3(0) = cxp 

Ax/6) - s(x/Q)e~ 2x ' 3 



G c {x) = 2- 
it is readily checked that 
/(0)=s(0/2)exp 



c(x/2) 



(5.23) 



dx G c (x) , 2 x 

T^r sm( ^ (i7r+0)) 

(5.24) 

minimally satisfies the above condition. To add the 
bound state structure, we note that the two fermions 
c, d form a bound state of mass \/3m. This implies a 
pole in the c-d S-matrix at = iu = in/3, and so a pole 
in f(0\ 2 ) at O2 — 81 + in/3. Hence we multiply f(8) by 
l/(c(0) - 1/2), leading us to the form for „f^(0 1 ,0 2 ): 

nfcd (0i, 02) = iA G (5 ac 6 bd - S ad 6 bc ) 

x ( e (f>i+02)/2 _ , 1 )M c -(9l+e 2 )/2- ) s (0i2/2) 

c(0i2) - 1/2 



x expl 



dx Gr[x) o / x 
~M W{ 2n 



sin 2 ( — (mt + 0i 2 ))], 



(5.25) 



where iA G is some normalization with mass dimension 

The phase of Aq is determined through the hermiticity 
condition 



m/c1(0i,02)* = i*,0i-in). 
This implies that Aq is real. 

Calculation of 



(5.26) 



Again we begin by identifying the group theoretical 
structure of 



^ /a/3 (0i, 02) 



(Ca ab ) a0 f tt (e 1 ,O 2 ), 



(5.27) 



where C is the charge conjugation matrix introduced in 
Section 4. (Ca' ab ) Q , ( g is not only anti-symmetric in a, b but 
also in a, (3 as it must be if the kinks are to couple to G° b 
fl36| |. Now (C<j ab ) al 3 is not the only available invariant 
tensor. One also has (a ab C) a /3, anti-symmetric combina- 
tions built up out of 7 a C7 b or C^C^C, or some com- 
bination of all three (but not a ab as this choice violates 
obvious U(l) conservation). Perhaps the most natural 

ab 

X0 



choice is to make ^rt explicitly C-symmetric: 



= ({Co ab ) aP + (a ab C) a p)f,(8 ll 6 2 ). (5.28) 

However this forces M / Q ^ a_1 ' )( ' 2a ' 1 to zero and this is not 

consistent with triality, i.e. all M / C ^ a 1 ^ 2a ^ do not van- 
ish. So we choose as above. Note that we could equally 
well have chosen (a ab C) a p. That our choice is not C- 
symmetric is not surprising as the theory is not triv- 
ially C-symmetric. With kink-fcrmion scattering, a C- 
transformation changes the sign of selected amplitudes 



(seepL20p. 

As before Lorentz invariance demands 



fu 



i) = (e 



(01+02)72 



(_l)M e -C*+fe)/a )/(6)l2) 



(5.29) 



Then by |5_2|, f(0 12 ) must satisfy 

f 36* in 



/(-0) = -u 2 (0) 



Ani -39 6 
= -(a 2 (e)-a 3 (e))f(0) 



-/(0), 



(5.30) 



where the last line follows from 4.26. That 5.30 has ex- 



actly the same form as 5.22 is a non-trivial consequence 



of triality. The periodicity axiom also takes the same 
form as before: 



f(-0) = f(0-2ni). 



(5.31) 



As the bound states of two fermions are identical to that 
of two same chirality kinks, the bound state structure of 



ftfaft is tne same as p. fed- Hence 



f ,f^(8i,0 2 )= i B(Ca ab ) af3 
,(0i+0 2 )/2 _ r_i\n e -(ei+6s 



x 



)/ 2) K012/2) 
'c(0 12 )-l/2 



x exp 



dxGJ 

~lU Smi ^ {in + ei2)) . 



(5.32) 
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What is left is to determine the relative value of B to Aq 
via triality. An overall normalization is not so interesting 
when computing two-point correlators as the physics so 
encoded is not universal. 

To determine the normalization we focus upon a par- 
ticular current component Gj 2 . Gj 2 is given by 



Gj 2 oc <9 x </>i 



(5.33) 



where <j)\ is one of the four Cartan bosons. Under the 
triality transformation of 2.25, Gq 2 is transformed into 



^12 . 1 (rili . n Zi , ^(56 | ^<78\ 

^ 2^ + o 



(5.34) 



Now focus upon the form factor where \ot) = 

(1/2, 1/2, 1/2, 1/2) and \0) = (-1/2, -1/2, -1/2, -1/2). 
Under a naive counting of quantum numbers, A a and Ap 
are transformed into Dirac fermions: 



1 



iAx) 



A a ^ -^A 2 
A p ^ —(A 2 -iA 1 ). 



Hence we have 



Jaf3 



1 rl2 
o J 12 ■ 



(5.35) 



(5.36) 



This fixes B = A G /2. 



E. One Particle Form Factors for SO (8) Currents 



Only the rank-2 tensorial states couple individually to 
the G° h 's. Denoting the Faddeev-Zamolodchikov opera- 



tors for the bound states by As cc n we have (as in 4.34) 

*fl&}4{cd}(0) = res 5=0 A c {0 + S + iir/6)A d (d - in/6), 

(5.37) 

where A c ,A d are Majorana fermions. 9^f^ , the ampli- 
tude to form the bound state At cd y from A c and Ad, is 

given by the residue in at iu]^: 



n {cd} n cd 
Med 9{cd} 



and so is readily computed to be 

: r(2/3K 1/2 



Jed} 



(2V37T 



r(i/6)' 



(5.38) 



(5.39) 



Note that the value of gj,^ given here is ambiguous up to 
a phase. However the phase can be pinned down through 
an appeal to hermiticity. 



Using 5.37, we can write down 



i*f{cd}(Q) = iA G (5 ac S bd - S ad 5 bc ) 



x(e 6 



x exp 



-1/2 



dx G c {x) 2 



o x s(x) 



(5.40) 



The overall phase (up to a sign) has been fixed by the 
constraint, M / { a c 6 d }* = -ju/fid}) arism S from hermiticity. 
Although the overall sign is uncertain, the computation 
of correlators does not depend on its value. 



F. Two Particle Form Factors for Gross-Neveu 
Fermion Operators 

At the two particle level, two kinks of opposite chirality 
will couple to a right/left fermion field, tjj±. This leads 
to form factors, ±f®g(6i, 6 2 ), defined by 

±/^(0i,0 2 ) = (^(0)^(02)^(01)). (5.41) 
Covariance under SO(8) demands ±f£p takes the form 

±f: P (0i,e 2 ) = [c 1 (c 1 a ) aP + ^0)^(9,, e 2 ), 

(5.42) 

where c\ and c 2 are constants. (That c\ and c 2 are 
not more generally independent functions of B\ and 9 2 
is easily seen; the constraints |5.2| and 5J3 do not al- 
low it.) c\ and c 2 can be fixed easily. Suppose \a) = 
(1/2,1/2,1/2,1/2) and \0) = (1/2,-1/2,-1/2,-1/2). 
Then by rewriting the Majorana fields as Dirac fields, 
one readily finds that 



±fa/3 — i±fa/3- 

This in turn forces c 2 = 0. We can then set c\ 
are uninterested in an overall normalization. 
Lorentz invariance demands 



f±(0 



1,02) 



D ±(01+02)/4 



/o(M- 



(5.43) 
1 as we 

(5.44) 



fo is then constrained through 5.2 



fo(-6) = -S 2 (in-6) 
= -S 2 (9)f (9), 



2m + W 
5ni - 30 



/o(0) 



where S 2 is given in 4.22, and through n 



T»(±%)/o(« - 27ri) = / o (-0), 



(5.45) 



(5.46) 



where ±R aa = ±i is the braiding phase. Given that 
S 2 (6) can be written as 
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S 2 (B) 



exp 



q X Ln 



G f (x) = 



e -7x/G + 2c(g/g) 

c(x/2) 

we see that /o takes the form 



/o(0) = exp 



dx Gf(x) . ? , x , . ... 
o a; s(x) Ztt 



(5.47) 



(5.48) 



As two kinks of opposite chirality form a fermionic bound 
state, i/ojg should have a pole at 2 = #i + iit®^ = 
6»i + i27r/3. Thus ±f% becomes, 



for f ±(9i,02) and f ±(61,62)- By triality we expect these 
forms to match that of ±/afl an d onr y f ±(61,62) does. 

We must also consider ±fg a , the form factor obtained 
from ±/„0 by reversing the order of particles. By trial- 
ity both ±fap an d ±fp a are expressible in terms of form 
factors of type ±f£p- Because is symmetric in a, (3, 
we must then have 



/a ffl 
af3 — ±Jl3a' 

and so ±f$ a (0i,02) = (Cr) a 0f±(0i,02)- 
Lorentz invariance demands that 



f±( 



(5.54) 



(5.55) 



±f 0(3(61,82 



) = A F (C 1 a ) a0 



e ±(e 1 +e 2 )/4 
c(0 12 ) + l/2' 



This ansatz is again constrained by 5.2 and 5.3 
junction with the kink-fermion S-matrix 1.20] : 



m con- 



exp 



dxGf(x) . 0, x ,. „ 
TV sm (T" ( m + e i*)) 

q x s(x) 2ir 



(5.49) 



where Ap is some normalization with mass dimension 



in 



1/2 



To determine the phase of Af we again employ a her- 
miticity condition: 



±/2/j(0i.ftO* = ±/!a(fc-«r,0i-t7r 



(5.50) 



This implies the phase of Ap is e ±47r / 4 . Scaling this phase 
out leaves us with the final form for ±f£p : 



/o(0-27ri) = -£ 2 (0)/o(0); 
fo(6-2m) = f o (-0). 



(5.56) 



That these equations match 5.44 and |5.45 is a reflection 
of triality. Given the bound structure of ±f"g 1S the same 
as ±/a/3i we can wr it e down 



= -Ape^^CYU 



e ±{0i+6 2 )/4 
C(012) + 1/2 



±fZp(0M = Ape^l^C^U^—. -^x 



0(6x2) + 1/2' 



x exp 



dx Gt(x) . o, x ,. n , , 
rV sin ^- ^ + ^2 

x s(x) 2w 



, (5.57) 



exp 



dx Gf(x) , o, x ,. n 

' 7Vsm 2 ( — (z^ + 12 )) 

q a; s(x) 27T 



where the normalization relative to ±/qa has been fixed 
(5.51) using triality. 



G. Two Particle Form Factors for Kink Operators 

At the two particle level, a right/left kink field, V±> 
will couple to a kink Ap(0) of chirality opposite to a to- 
gether with a Gross-Neveu fermion. So we consider form 
factors given by 

±WM = (r ± (Q)Ap(02)A a (9 l )). (5.52) 

For ±f™g to be invariant under 50(8), it must take the 
form 

±WM = (c 1 a ) aP f ± (0 1 ,02) + (i a c) aP ~f ± (0 1 ,02). 

(5.53) 

We can set f ±(61,62) to 0. The constraints upon the 
form factors 5.2 and 5.3 lead to different functional forms 



H. One Particle Form Factor for Kink Fields 

The kink field, will couple to the anti-kink C a pAp. 
Hence we consider the form factor, 



±f$(6) = (^(0)4,(0)), 
which must take the form 

±fp(6) = c±e ±e ' 2 C aP . 



(5.58) 



(5.59) 



To determine the constant, we use the two particle form 
factor ±f a afj . Let \a) = (1/2,1/2,1/2,1/2) and \(3) = 
(-1/2,-1/2,-1/2,-1/2). Ap can be written as 

ig%A (6) = res s =oA a (6 + S + tu a ^)A 7 (6 - m^), 

(5.60) 
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where a = 1, (7) = (1/2,-1/2,-1/2,-1/2). Again the 
it's mark out poles in S-matrix indicative of bound states. 
Here they are given by 



U a 7, = U 1 a = 7T — U 1 , — 7r/3. 

7/3 a/3 a(3 1 



(5.61) 



g@ can be determined up to a phase from the fermion- 
kink S-matrix as before to be 



9 ^ = 2V WW) 



1/2 



(5.62) 



Given that, 



i$£r±/|(0) = res o=0 ± f« a (0 - in/3, 9 + 5 + in/3), 



(5.63) 



together with the hermiticity constraint, ±/™(#)* 
(+i)±f£{0), we find (up to a sign) 



= iA G (6 ac 5 bd - 5 a d6bc)fn(0i, O2), 
,C(9 1 ,e 2 )^(Gf(Q)A p (e 2 )A a (6 1 )) 

where 

MOiM = ( e (^+^)/ 2 - (_i)M e -(»i+« a )/a) 
*(»ia/2) 



(5.67) 



c(M - 1/2 
x exp 



da; G c (x) . 2 . x 



G c (x) = 2 



q X s(x) 27T 

c(x/6) - s(x/6)e" 2x / 3 
c(a;/2) 



sin 2 ( — (t7r + la )) 



-1/2 



x exp 



^!i^ s2(x / 6 ) 

x s(x) 



(5.64) 



I. One Particle Form Factor for the Gross-Neveu 
Fermions 

At the one particle level, the fermion field, ip± will 
couple only to A a {9). Hence we have 

± fS = (V4(0)A(^)> = d ± e ±e / 2 S ab . (5.65) 

To determine the constant d±, we use triality. Fixing 
a = 2, it is easy to show under the triality transforma- 
tion |!2f , 



±fa - ^(i/a + ±fS)> 

where a = (1/2, 1/2, 1/2, 1/2) and so d± = c±. 



(5.66) 



(5.68) 

For the Gross-Neveu fermions, ^>±, and the kinks, ip±, 
we have 

= A F e ± ™'\C 1 a ) a0 f±{9 1 ,9 2 ), 
= -A F e ± ™/\C 1 a ) aP f ± (e 1 ,6 2 ), 



where 



x exp 



Gf(x) 



e ±(e 1 +e 2 )/4 

c(9 12 ) + 1/2 

dxGf(x) . 2/ x 

2c(x/6) + e~ 7x ' 6 
c{x/2) ' 



(5.69) 



(5.70) 



J. Summary 

Here we summarize the results of this section for quick 
reference: 

Two Particle Form Factors: 

For the 5*0(8) currents, G° b , we have 



l }^{e 1 ,e 2 ) = {Gf{Q)A h {e 2 )A a {e x )) 



One Particle Form Factors: 



For the currents, G^ b , we have 



,ff cd} (6) ee (G ab (O)A {cd} (0)) 

= iA G {S ac 5 bd - 6 ad 6 bc )(e e - (-l) M e" 
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x exp 



dx G c ) > 



lo x s(x) 
and for the fermions, and kinks, ip^_, 

±/ b a = (r ± (0)A b (9)) =c±e ±e ' 2 8 ab , 
±/| = (^(Q)A a (9)) = c±e ±e / 2 C aP , 
where the constant c± is 

-1/2 



(5.71) 



(5.72) 



V3 



x exp 



r(7/6)y 

o a; s(x) 



(5.73) 
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